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Abstract. Generally the gas metallicity in distant galaxies can only be inferred by using a few prominent emission 
lines. Various theoretical models have been used to predict the relationship between emission line fluxes and 
metallicity, suggesting that some line ratios can be used as diagnostics of the gas metallicity in galaxies. However, 
accurate empirical calibrations of these emission line flux ratios from real galaxy spectra spanning a wide metallic- 
ity range are still lacking. In this paper we provide such empirical calibrations by using the combination of two sets 
of spectroscopic data: one consisting of low-metallicity galaxies with a measurement of [Oni]A4363 taken from the 
literature, including spectra from the Sloan Digital Sky Survey (SDSS), and the other one consisting of galaxies in 
the SDSS database whose gas metallicity has been determined from various strong emission lines in their spectra. 
This combined data set constitutes the largest sample of galaxies with information on the gas metallicity available 
so far and spanning the widest metallicity range. By using these data we obtain accurate empirical relations be- 
tween gas metallicity and several emission line diagnostics, including the R23 parameter, the [Nn]A6584/Ha and 
[Oiii]A5007/[Nii]A6584 ratios. Our empirical diagrams show that the line ratio [Om] A5007/[On]A3727 is a useful 
tool to break the degeneracy in the R23 parameter when no information on the [Nn]A6584 line is available. The line 
ratio [Nem] A3869/[On]A3727 also results to be a useful metallicity indicator for high-z galaxies, especially when 
the R23 parameter or other diagnostics involving [Oni]A5007 or [Nll]A6584 are not available. Additional, useful 
diagnostics newly proposed in this paper are the line ratios of (Ha+[Nll]AA6548,6584)/[Sll]A6720, [Oni] A5007/H/3, 
and [On]A3727/H/3. Finally, we compare these empirical relations with photoionization models. We find that the 
empirical i?23-metallicity sequence is strongly discrepant with respect to the trend expected by models with con- 
stant ionization parameter. Such a discrepancy is also found for other line ratios. These discrepancies provide 
evidence for a strong metallicity dependence of the average ionization parameter in galaxies. In particular, we find 
that the average ionization parameter in galaxies increases by ~ 0.7 dex as the metallicity decreases from 2 Zq 
to 0.05 Zg, with a small dispersion. This result should warn about the use of theoretical models with constant 
ionization parameter to infer metallicities from observed line ratios. 
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1. Introduction 

The gas metallicity is one of the most important tools 
to investigate the evolutionary history of galaxies. This 
is because the gas metallicity of galaxies is basically de- 
termined by their star-formation history. Recent observa- 
tional studies allowed the investigation of the gas metal- 
licity even in high-z galaxies beyond z = 1, such as 
Lyman-break galaxies (e.g., Teplitz et al. 2000a, 2000b; 
Pettini et al. 2001), submillimeter-selected high-z galax- 
ies (Swinbank et al. 2004), and so on (see also, e.g., 
Kobulnicky & Kewley 2004; Savaglio et al. 2005; Maier 
et al. 2006; Liang et al. 2006; Erb et al. 2006). Such obser- 
vational insights on the metallicity evolution of galaxies 
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are now giving constraints on the theoretical understand- 
ings of the formation and the evolution of galaxies (e.g., 
Bicker et al. 2004). 

However, metallicity measurements in distant galaxies 
are not straightforward. Information on the gas temper- 
ature is required for a precise determination of the gas 
metallicity, but the gas temperature can be accurately in- 
ferred only when the fluxes of auroral emission lines such 
as [Oiii]A4363 and [Nn]A5755 are known, and these are 
generally too weak to be measured in faint distant galax- 
ies. The measurement of the auroral emission lines is diffi- 
cult even for galaxies in the local universe especially when 
the gas metallicity is high, because the collisional excita- 
tion of the auroral transitions is suppressed due to efficient 
cooling through far-infrared fine-structure emission lines 
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(e.g., Ferland et al. 1984; Nagao et al. 2006a). Therefore, 
in general we have to rely on some relations between gas 
metallicity and flux ratios of strong emission-lines to esti- 
mate the gas metallicity in most galaxies. Extensive stud- 
ies have been performed to calibrate such metallicity di- 
agnostics by using only strong emission lines. One of the 
most frequently used metallicity diagnostics is the R23 pa- 
rameter, defined as 



-R23 = 



J 1 ([OII]A3727)+F([OIII]A4959)+F([OIII]A5007) 
F(H/3A4861) : 



(1) 



where F([On] A3727), F([Om]A4959) and so on denote the 
emission-line fluxes of [On]A3727, [Oiii]A4959 and so on, 
respectively. The R23 was proposed by Pagel et al. (1979), 
and its calibration to the oxygen abundance has been 
improved by various photoionization model calculations 
(e.g., McGaugh 1991; Kewley & Dopita 2002). 

One serious problem of this indicator is that a certain 
value of i?23 has two different solutions, a low-metallicity 
solution and a high metallicity one. Therefore additional, 
or alternative, diagnostics aimed at removing the R23 de- 
generacy have been proposed (e.g., Alloin et al. 1979; 
Denicolo et al. 2002; Kewley & Dopita 2002; Pcttini & 
Pagel 2004). However, most of these methods exploit the 
[Nn]A6584 line, which has the problem of being very weak 
at sub-solar metallicities (hence difficult to measure) and 
the problem of being rapidly shifted outside the spectral 
band of many surveys at high redshift (e.g., unusable be- 
yond z <~ 0.5 in optical spectra). On the theoretical side 
various models have been presented, which provide the ra- 
tios among the most prominent emission lines as a function 
of metallicity (e.g. Kewley and Dopita 2002). However, 
model predictions strongly depend on the assumed phys- 
ical parameters of the ionized gas, and in particular on 
the ionization parameter [U = Qn/(cnn), where <E>h is the 
surface flux of hydrogen-ionizing photon and nn is the gas 
density]. As a consequence, an accurate correspondence 
between individual diagnostics (line ratios) and metallic- 
ity cannot be established, because of the lack of informa- 
tion on the physical conditions of the gas. Summarizing, 
many gas metallicity diagnostics proposed so far are ei- 
ther ambiguous or unusable when applied to the spectra 
of distant galaxies. 

The goal of this paper is to obtain accurate, empirical 
calibrations between metallicity and individual diagnos- 
tics involving a few strong emission lines, which can be 
applied to the spectra of distant galaxies. In particular, 
we re-calibrate diagnostics already proposed in the past, 
but we also propose new diagnostics which appear partic- 
ularly suited for distant galaxies. This work is obtained by 
combining two large data sets. The first one is composed of 
recent spectroscopic observations of low-metallicity galax- 
ies [7.0^12+log(O/H)^8.5], whose metallicity is accu- 
rately determined through the [Oiii]A4363 line (§§2.1). 
This dataset consists of two subsamples; one is taken from 
the database of the Sloan Digital Sky Survey (SDSS; York 
et al. 2000; Strauss et al. 2002) (§§2.1.1) and the other is 
taken from the literature (§§2.1.2). The second data set 



is a subsample of galaxies in the SDSS database, whose 
metallicity has been derived by Tremonti et al. (2004) 
[8.2^12+log(0/H)^9.2; see §§2.2]. These combined data 
sets provide the largest sample of galaxies with informa- 
tion on the gas metallicities and spanning more than 2 dex 
in metallicity. 



2. Data 

2.1. Spectroscopic data of low-metallicity galaxies with 
a [OIIIJX4363 measurement 

2.1.1. SDSS data (sample A) 

The gas-phase oxygen abundance is well determined when 
the flux of [Om]A4363 is measured (e.g., Osterbrock 1989). 
Although such measurements have been performed for 
more than a hundred low-metallicity galaxies, simple com- 
pilation of those earlier results may introduce some un- 
expected biases and uncertainties. This is because the 
data were collected by various (heterogeneous) observa- 
tions with different properties (aperture size, wavelength 
resolution, and so on) and because the method of calculat- 
ing the oxygen abundance is different for different authors. 
Recently, Izotov et al. (2006b) reported their systematic 
measurements of the oxygen abundance for low-metallicity 
galaxies in the SDSS Data Release 3 (DR3; Abazajian 
et al. 2005) by using the [Oiii]A4363 emission-line flux. 
The extinction-corrected emission-line fluxes of galaxies 
with a measurement of the oxygen abundance provided 
by Izotov et al. (2006b) are the ideal data for the em- 
pirical calibration of metallicity diagnostics, because the 
data were obtained and measured in a homogeneous way 
and because the oxygen abundance is also calculated with 
a common method. The number of spectra analyzed by 
Izotov et al. (2006b) is 309. Among them, we use the 
data with a relatively small error in the oxygen abun- 
dance [A(log§) < 0.05] (146 spectra). Here we adopt the 
uncertainty [A(log§)] given in Table 2 of Izotov et al. 
(2006). Since some spectra in the database of Izotov et al. 
(2006b) are duplicated for the same objects, the number 
of galaxies with A(log§) < 0.05 is 139. Hereafter we call 
this sample "sample A" . 

However, this sample has two problems when used to 
accurately calibrate metallicity diagnostics. First, most 
of galaxies in sample A have a relatively high oxygen 
abundance, and only 6 of them have 12+log(0/H) < 7.6. 
Therefore the statistical reliability of the empirical cali- 
bration of metallicity diagnostics would be extremely poor 
at 12+log(0/H) < 7.6 if using only this sample. Second, 
in sample A, there is the remarkable tendency for lower- 
metallicity galaxies to have lower redshift. In Figure 1, 
the oxygen abundance of galaxies in sample A is shown 
as a function of redshift. The origin of this apparent ten- 
dency is likely due to the fact that SDSS is not a volume- 
limited survey; that is, galaxies at higher redshift have 
preferentially higher luminosity, and thus higher metallic- 
ity. In particular, all of the galaxies with 12+log(0/H) 
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< 7.6 are at z < 0.02. This means that the [On]A3727 
flux cannot be measured for the latter galaxies due to 
the limited wavelength coverage of the SDSS spectroscopy 
(A bs *S 3800A). Therefore, if using only sample A, we 
could not calibrate the diagnostics involving [On]A3727 
[R 23 and F([Nll]A6584)/F([Oll]A3727)] in the metallicity 
range 12+log(0/H) < 7.6. This is a serious problem, be- 
cause i?23 is one of the most frequently used metallicity di- 
agnostics and thus should be calibrated in a wide metallic- 
ity range. In conclusion, the accurate calibration of various 
metallicity diagnostics in a wide metallicity range cannot 
be achieved by using only sample A. We therefore collected 
additional data of [Oin] A4363-detected galaxies, which are 
described in the following subsection. 

2.1.2. Other data from literature (sample B) 

In order to increase the number of low- metallicity galaxies 
with a measurement of the oxygen abundance, we com- 
piled the reddening-corrected emission-line flux data of 
galaxies with a [Oiii]A4363 measurement from the litera- 
ture. The sample of compiled galaxies is given in Table 
1. For the objects whose spectroscopic properties have 
been reported by more than one paper independently, we 
chose the one with higher signal-to-noisc ratio. When both 
the spectroscopic properties of the whole galaxy and of 
parts of it have been reported, we compiled both of them 
(e.g., Mrk 116). Consequently, the number of the com- 
piled objects is 157. To minimize possible systematic er- 
rors owing to the different methods on the calculation of 
the oxygen abundance, we re-calculate their oxygen abun- 
dance by adopting the same method used for sample A 
(Izotov et al. 2006b). The re-calculated R23 parameter, 
gas density of [Sn] -emitting region [nu(S + )], gas temper- 
ature of [Om]-emitting region [i(0 2+ )] and oxygen abun- 
dance [12+log(0/H)] are given in Table 1, along with the 
reference to the data of the emission-line flux ratios. To 
calculate R23, we did not use F([Oiii]A4959) but calcu- 
late the ratio of [F([On]A3727) + 1.327xF([Om]A5007)] 
/F(H/3) since only the flux of [Om]A5007 (without that of 
[Oiii]A4959) is given in some reference papers. 

To check whether our adopted method causes possible 
systematic difference in the oxygen abundance from the 
values given in the original references, we compare the 
oxygen abundances re-calculated by us and those given in 
the original papers in Figure 3. Apparently, there is no 
systematic difference between our results and the results 
given in the literature. The mean and the RMS of the dif- 
ference, [12+l0g(0/H)this work] - [12+log(0/H)i iterature ], 
are +0.001 and 0.041, respectively. This mean value of 
the difference is smaller than the typical error of the 
re-calculated oxygen abundance. Among the 157 objects 
given in Table 1 and plotted in Figure 3, we use only 
data with A(log§) < 0.05. This constraint results in a 
sample of 120 objects, that is hereafter called "sample 
B". Note that the mean and the RMS of the difference, 
[12+log(0/H)this work] - [12+log(0/H) lit oraturo], for sam- 



ple B are -0.008 and 0.038. Again the mean value of the 
difference is smaller than the typical error on the oxygen 
abundance. 

2.2. Spectroscopic data of high-metallicity galaxies 
(sample C) 

For high-metallicity galaxies, we referred to the oxygen 
abundance derived by Tremonti et al. (2004), who de- 
rived the metallicities of ^53000 galaxies in the SDSS 
database. They used the fluxes of many strong emission 
lines ([On]A3727, H/3A4861, [Om]A5007, Ha, [Nn]A6584, 
[Sn]A6717, and [Sn]A6731) and comparing them with pho- 
toionization models (Ferland et al. 1998). Although they 
presented the results of their analysis on the spectra of 
the SDSS Data Release 2 (DR2; Abazajian et al. 2004), 
they also provide the results of their recent analysis on the 
spectra from Data Release 4 (DR4; Adelman-McCarthy 
et al. 2006) on their web site 1 . Their estimate of the oxy- 
gen abundance does not rely only on a single metallic- 
ity diagnostic flux ratio, but uses all the optical promi- 
nent emission lines (see also, Chariot & Longhetti 2001; 
Brinchmann et al. 2004). Therefore, among galaxies with- 
out [Om]A4363 flux, their sample is currently the best one 
in terms of both sample size and reliability. 

The oxygen-abundance catalog of the SDSS DR4 
galaxies contains 567486 objects. The objects in the cata- 
log are classified into the five classes; "star-forming galax- 
ies", "low S/N star-forming galaxies", "composite", "ac- 
tive galactic nuclei (AGNs)" , and "unclassificable" . We re- 
ferred only to the galaxies belonging the first class (141317 
objects). The emission- line fluxes of these galaxies are ob- 
tained from the emission-line flux catalog of the SDSS 
DR4 galaxies provided on the same web site as the oxygen- 
abundance catalog. The emission-line fluxes given in this 
catalog were measured from the stellar-continuum sub- 
tracted spectra with the latest high spectral resolution 
population synthesis models by Bruzual & Chariot (2003), 
and thus more reliable than the flux data provided on the 
SDSS Data Archive Server. Since the emission-line fluxes 
given in their catalog are not corrected for dust extinc- 
tion, we corrected them by using the Balmer decrement 
method with the reddening curve of Cardelli et al. (1989). 
We then removed the duplicated objects and the objects 
observed in some problematic plates (see the SDSS web 
page 2 ) from the cross sample of the oxygen-abundance 
catalog and the emission-line flux catalog. Then we select 
only objects satisfying all of the following five criteria: 

1. The rcdshift is higher than 0.028. 

2. Both Ha and H/3 emission lines have S/N > 10. 

o , F(fOIIIlA5007) 0-61 , i o 

6 - 10 § F(H/3) > log F(|Nii^6584) ~ + 

4. The fiber aperture covers at least 20% of the total g'- 
band photons. 

1 http:/ /www. mpa-garching.mpg.de/SDSS/ 

2 http://www.sdss.org/dr4/ 
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5. The uncertainty on the estimated stellar mass is less 
than 0.5 dex (i.e., logM 9 7. 5 - logM 2 .5 < 0.5, where 
M97 5 and M2.5 are the 97.5th and 2.5th percentile of 
probability distribution of the estimated stellar mass; 
see Kauffmann et al. 2003b for more details). 

The first criterion ensures the accurate measurement of 
.F([Oii] A3727). The second criterion selects emission-line 
galaxies. The S/N values are taken from the emission-line 
flux catalog; however, the errors in this catalog are under- 
estimated with a factor of ^2 (see the web site for more 
details). The third criterion is required to reject AGNs 
from the sample, following Kauffmann et al. (2003a), al- 
though the AGN removal had been already examined in 
the process of checking the classification in the oxygen- 
abundance catalog mentioned above. The fourth criterion 
is required to avoid significant aperture effects on the flux 
ratios. Kewley et al. (2005) reported that the introduced 
systematic error of the metallicity determination reaches 
up to ^40% when the fiber covers only 20% of the to- 
tal <7'-band photons (see also Trcmonti et al. 2004). The 
fifth criterion is not directly relevant to this study, but is 
required in our companion paper (Nagao et al. 2006b). 

Note that in sample C we do not put constraints on 
A(log§). This is to prevent sample C to be devoid of 
galaxies at metallicities 12+log(0/H) < 7.6 so that sample 
C and sample A+B have some overlap in terms of metal- 
licities. Indeed, due to the luminosity-metallicity relation, 
low metallicity galaxies in sample C are on average fainter 
and therefore tend to have larger errors. However, the av- 
eraged gas metallicities in sample C are most likely reli- 
able even at low metallicities, thanks to the large number 
of objects in the sample. Our final sample of emission-line 
galaxies consists of 48497 objects, which is hereafter called 
"sample C" . The redshift distribution of galaxies in sam- 
ple C is shown in Figure 4. Its median value is 0.085, while 
the mean and the RMS are 0.092 and 0.040, respectively. 
The distribution of the oxygen abundance of galaxies in 
sample C is shown in Figure 5. Its median value is 9.016, 
while the mean and the RMS are 8.976 and 0.166. The 
means and the RMSs of the oxygen abundance of samples 
A, B, and C are summarized in Table 2. 

3. Results 

In Figures 6 and 7, we plot emission-line flux ratios 
for the galaxies in samples A, B, and C. To avoid 
noisy objects in sample C, we consider only those with 
S/N > 10 (cataloged value) for all the related emission 
lines (e.g., H/3, [On]A3726, [On]A3729 and [Om]A5007 
for the case of #23)- In addition to R 2 ^ 1 all the 
other flux ratios investigated here are metallicity-sensitive 
flux ratios and sometimes regarded as metallicity di- 
agnostics (see, e.g., Kewley & Dopita 2002; Pettini & 
Pagel 2004; Kobulnicky & Kewley 2004). Among them, 
F([Om]A5007)/F([Oll]A3727) is sensitive also to the ion- 
ization parameter and thus it has not been regarded as 
a good metallicity diagnostic flux ratio (see Kewley & 



Dopita 2002). Instead, this flux ratio has been used to 
investigate the ionization parameter, and has been some- 
times used in the following form: 

_ F([OIII]A4959)+F([OIII]A5007) 
32 ~ F([OII]A3727) ( ' 

(e.g., Kobulnicky & Kewley 2004). 

The data sequences in the diagnostics-metallicity dia- 
grams are mostly continuous for different samples, and ac- 
cordingly the whole sample shows clear relations between 
various metallicity diagnostics and the oxygen abundance. 
Since there are no apparent systematic differences in the 
diagnostics-metallicity sequences between sample A and 
sample B, we combined these two samples and identi- 
fied as "sample A+B" hereafter in order to improve the 
statistics at low metallicities. The statistical properties 
of the sample A+B are given in Table 2. The diagram 
of F([Nn]A6584)/F([Sll]A6720) versus the oxygen abun- 
dance shows an apparent discontinuity between sample 
A+B and sample C, where F([Sn]A6720) denotes the sum 
of F([Sll]A6717) and F([Sll]A6731). We will discuss the 
issue of this discontinuity in §§4.1. 

The diagram of i? 2 3 versus the oxygen abun- 
dance shows a H-shaped distribution with a peak at 
12+log(0/H) ~ 8.0. This is consistent with the previous 
studies on the empirical relation between i? 2 3 and the oxy- 
gen abundance based on smaller samples of observational 
data (e.g., Edmunds & Pagel 1984; McGaugh 1991; Miller 
& Hodge 1996; Castellanos et al. 2002; Lee et al. 2003a; 
Bresolin et al. 2004, 2005; Pilyugin & Thuan 2005). As 
discussed in §§4.3, however, this appears to be systemati- 
cally different from previous predictions of photoionization 
models. 

To investigate the relation between the flux ratios and 
the oxygen abundances quantitatively, we calculate the 
means and the RMSs of the flux ratios of galaxies within 
bins of oxygen abundance. For sample A+B, we calculate 
them in the range 7.05 < 12+log(0/H) < 8.55 with a bin 
width of A[log(0/H)] =0.1, except at lowest and high- 
est oxygen abundances, where the bin width is wider (i.e., 
7.05 < 12+log(0/H) < 7.45 and 8.35 < 12+log(0/H) 
< 8.55) due to the small number of sources in these ranges. 
All of the metallicity bins contain at least 6 galaxies. 
The results are given in Table 3. We also calculate the 
mean and the RMS of flux ratios for sample C in the 
range 8.15 < 12+log(0/H) < 9.25 with a bin width of 
A[log(0/H)] = 0.1 dex. The results are given in Table 4. 

The calculated mean and the RMS of the flux ratios 
for each metallicity bin are shown in Figures 8 and 9. We 
then fit the observed sequences between flux ratios and 
oxygen abundance with polynomial functions in the range 
7.05 < 12+log(0/H) < 9.25 (or 0.02 < Z gas /Z Q < 4), 
and the results of the fits are also shown in Figures 8 
and 9. We decided to fit 3rd-order polynomial functions 
for the binned data, not for the individual data, in order 
to avoid giving too much weights to the high metallicity 
range (where most of the data are). In Table 5, the coef- 
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ficicnts of the best-fit polynomial functions are provided, 
according to the formula 

log (Fx/F 2 ) = Yl a N [log(Z/Z e )} N (3) 

where Fi / F 2 is the line flux ratio (or the i? 2 3 parameter) 
and N = (0, 1, 2, 3). Here we adopt 12+log(O/H) = 8.69 
for Zq (Allende Prieto et al. 2001). For the convenience of 
the reader, in Table 6 we also give the coefficients of the 
best-fit polynomial functions in the following form: 

log (F 1 /F 2 ) = ]T M12 + log(0/H)] w . (4) 

JV 

The expected uncertainties on the derived metallicities 
from the diagnostic flux ratios calibrated here can be 
estimated using the RMS values of the flux ratios given 
for each mass bin (Tables 3 and 4). By looking at the 
RMS plotted in Figures 8 and 9, we can recognize that, 
for instance, the diagnostic flux ratios give highly uncer- 
tain metallicities when F([Nii]A6584)/F([Oii]A3727) 
< 0.05, F([Oiii]A5007)/F([Oii]A3727) > 2, and 
J F([Nll]A6584)/ J F([Sll]A6720) £ 0.3. On the con- 
trary, the metallicity is well determined (AZ ^ 0.2 
dex) when F([Oiii]A5007)/F([Nii]A6584) £ 10 and 
F([Nii]A6584)/F([Oii]A3727) £ 0.05. 

Finally we recall that these relations are valid only in 
the range 7.05 < 12+log(0/H) < 9.25 (which is however 
much wider than in any previous work). We warn on the 
use of these relations outside such metallicity range, since 
it would rely only on their extrapolation. 

4. Discussion 

4.1. Consistency between sample A+B and sample C 

Before interpreting the results, we discuss on the 
consistency of the two main samples, i.e., galax- 
ies with (A+B) and without (C) [Oiii]A4363 mea- 
surements. As mentioned in §3, the relation between 
some emission-line flux ratios and the oxygen abun- 
dance is not smoothly connected between the two sam- 
ples (A+B and C), and this is especially significant 
for the flux ratios of F([Om] A5007)/F([Oll] A3727) and 
J F([Nll]A6584)/ J F([Sll]A6720), but is also seen in other 
cases (Figures 8 and 9). One of the possible reasons for 
this discrepancy is a systematic error in the estimate of the 
oxygen abundance for one (or both) of the two different 
methods, which in one case consists in using the gas tem- 
perature inferred through [Oiii]A4363 emission (see §§2.1) 
and in the other case is using all of optical strong emission 
lines (Trcmonti et al. 2004). 

Kobulnicky et al. (1999) investigated a possible sys- 
tematic error in the former method, that is, the gas tem- 
perature may be overestimated through the [Oiii]A4363 
emission and thus the oxygen abundance may tend to be 
underestimated accordingly. This is because the strength 
of the [Oiii]A4363 emission significantly depends on the 



gas temperature and thus spectra obtained by a global 
aperture toward a galaxy are biased towards higher 
gas-temperature Hn regions (see also Peimbert 1967). 
According to their analysis, the overestimation of the 
gas temperature could be more serious in low-metallicity 
systems and could reach up to AT C = 1000 — 3000K, 
which results in the systematic underestimation of the 
oxygen abundance of 0.05 - 0.2 dex. However, although 
this effect may partly account for the discrepancy of the 
metallicity dependence of F([Nii]A6584)/F([Sii]A6720), 
it goes in the opposite direction to account for the 
discrepancy seen in F([Olll]A5007)/F([Nll]A6584) and 
F([Om]A5007)/ J F 1 ([Oii]A3727). Therefore, the effect of the 
biased temperature measurement is not the dominant ori- 
gin of the discontinuities seen in Figures 8 and 9. 

A systematic error in the oxygen abundance may ex- 
ist in the method of Tremonti et al. (2004). They esti- 
mated the oxygen abundance by comparing photoioniza- 
tion models with some optical emission-line fluxes, which 
were measured on the spectra after subtraction of the 
stellar component. Although their method of the stellar- 
component subtraction is a sophisticated one and uses 
the most recent population synthesis models of Bruzual & 
Chariot (2003), it is not clear whether the measurement 
of emission lines lying on the deep and complex stellar 
absorption features is completely free from some possible 
systematic errors. A possible improper subtraction of the 
stellar absorption features may lead to systematic errors 
on the fluxes of Balmcr lines, which might result in a sys- 
tematic error in the estimation of the gas metallicity of 
galaxies in sample C. The subtraction of stellar absorp- 
tion features may be inaccurate also in sample A+B. For 
instance, in some earlier works the stellar subtraction was 
performed by simply assuming EW(H(3) a bs = 1A. This 
over-simplified assumption may introduce systematic er- 
rors in the derived gas metallicity and the emission-line 
flux ratios given in Table 1. Another possible source of un- 
certainty in the method of Trcmonti et al. (2004) is the use 
of the [Nn]A6584 flux and its comparison to models. Most 
photoionization models assume that the relative nitrogen 
abundance scales with the metallicity linearly when the 
primary nitrogen creation dominates, and scales quadrat- 
ically when the secondary nitrogen creation is dominant. 
However, the transition metallicity between the two modes 
is uncertain. An inaccurate value of the transition metal- 
licity (which is indeed uncertain) may lead to systematic 
errors in the estimation of metallicity especially at low 
metallicities, which could be one of the possible origin of 
the discrepancy seen in Figures 8 and 9. 

The discrepancy in the metallicity dependences of 
emission-line flux ratios may also be a consequence of 
the selection of spectroscopic targets. While galaxies in 
sample C are basically selected in terms of their appar- 
ent magnitude and thus not largely biased toward any 
specific population, galaxies in sample B could be bi- 
ased toward very strong emission-line galaxies (galaxies 
in sample A are in a composite situation; see Izotov et 
al. 2006). This is because the motivation behind most 



6 



Nagao et al.: Gas metallicity diagnostics 



of the original observations, such as the studies on the 
primordial helium abundance (see the original references 
given in Table 1), required very accurate measurements 
of emission-line flux ratios. For a given metallicity, galax- 
ies with stronger emission lines tend to be characterized 
by a higher ionization parameter, which may result into 
larger flux ratios of F([Olll]A5007)/F([Nll]A6584) and 
F([Om]A5007)/F([Oii]A3727), although the difference in 
the ionization parameter should not cause a significant 
difference in the ratio of F([Nll]A6584)/F([Sn]A6720). We 
will discuss the effect of the ionization parameter on the 
discrepancy further in §§4.3. 

Actually some or all of the above matters could con- 
tribute to the discontinuity in the metallicity dependences 
of emission- line flux ratios, and their discrimination or 
their accurate correction are not feasible. We thus sim- 
ply adopt the results of the fit described in §3 and not 
take the effects of the possible systematic errors into ac- 
count in the following discussion. However, it should be 
noted that this rather complex situation is caused by re- 
lying on two independent methods to measure the oxygen 
abundance. This problem will be solved if a large sample 
of galaxies with a wide range of the oxygen abundance 
is investigated by using a unique method throughout the 
concerned metallicity range. 

4.2. Comparison with previous empirical calibrations 

We compare the results of our calibrations with previ- 
ous empirical calibrations. In particular, in Figure 10, we 
compare the empirical calibrations of R23 derived by us 
with those obtained by Tremonti et al. (2004), Edmunds 
& Pagel (1984), and Zaritsky et al. (1994). While there is 
a reasonable agreement between our result and the result 
from previous calibration for the lower branch (Edmunds 
& Pagel 1984), there are some systematic discrepancies for 
the upper branch. We should in particular discuss the dif- 
ference between our calibration and that of Tremonti et al. 
(2004), since our calibration in the high- metallicity range 
is based on the metallicity of the SDSS galaxies (in sample 
C) derived by Tremonti et al. (2004). The calibration by 
Tremonti et al. (2004), which is provided only for the up- 
per branch, is clearly flatter than ours. This discrepancy 
may be ascribed to the combination of various possible 
factors. Our calibration also includes the fit of the new 
sample of [Oin] A4363-detected galaxies, which are not in- 
cluded in Tremonti et al. (2004), and this is certainly one 
of the reasons for the discrepancy. However, the latter is- 
sue cannot completely account for the discrepancy, since 
the Tremonti et al. (2004) calibration fails to reproduce 
the SDSS data at 12+log(0/H) < 8.5 (as shown in Figure 
10). It is likely that an additional source of the discrepancy 
is the different strategy of fitting the analytical function 
to the data. While we fit the third polynomial function 
to the binned data, Tremonti et al. (2004) fit the function 
to the whole sample of individual SDSS galaxies. Since 
the number of high metallicity galaxies [12+log(0/H) > 



8.5] is much larger than the low metallicity sub-sample 
[12+log(0/H) < 8.5] as shown in Figure 5, the analytical 
fit of Tremonti et al. (2004) is dominated by the high- 
metallicity part of the R23 diagram. Finally, the modest 
discrepancy at high metallicities [12+log(0/H) > 9] may 
be partly attributed to the difference in the sample se- 
lection criteria. As described in §3, we select the galax- 
ies in the sample C with S/N > 10 for all of the lines 
H/3, [On]A3726, [On]A3729, and [Om]A5007 (note that the 
[On] doublet lines are measured separately in the original 
catalog we used), while Tremonti et al. (2004) adopted 
the S/N criteria only for H/3, Ha and [Nn]A6584, not for 
[On]A3726, [On]A3729, and [Om]A5007. Our selection cri- 
teria may preferentially reject objects at high metallicities 
with respect to those of Tremonti et al. (2004), because 
forbidden lines such as [On] and [Oin] become weak when 
gas metallicity is high due to the suppressed collisional 
excitation mechanism (e.g., Ferland et al. 1984; Nagao et 
al. 2006a). This effect may result in our selective choice 
of objects with strong [On] and [Om] emission in a given 
metallicity bin, which could make our calibration to be 
steeper at high metallicities. Since the difference in the 
calibration between ours and that of Tremonti et al. (2004) 
is significant at I2+log(0/H) > 9, it is suggested that our 
calibration for the R23 may overestimate the gas metallic- 
ity at 12+log(0/H) > 9 by a factor of AZ ~ 0.1 dcx at 
12+log(0/H) ~ 9.1. 

The calibration of the diagnostic flux ratio 
F([Nn]A6584)/F(Ha) is especially important, be- 
cause wavelength separation of the two lines is small 
(i.e., not sensitive to dust reddening and requiring only 
small wavelength coverage) and thus it is used as a 
diagnostic of the gas metallicity of galaxies at z <> 2 
(e.g., Erb et al. 2006). In Figure 11, we compare the 
empirical calibrations of F([Nn]A6584)/f 1 (Ho:) derived 
by us, with those derived by Pettini & Pagel (2004) and 
Denicolo et al. (2002). Our result agree reasonably well 
with Denicolo et al. (2002) only at sub solar metallicity, 
and there is a systematic difference in the slope between 
our result and the result reported by Pettini & Pagel 
(2004). The latter discrepancy may be due to the reduced 
metallicity range of the sample of Pettini & Pagel (2004), 
indeed most of their objects are distributed within 7.7 < 
f2+log(0/H) < 8.5. However, the difference is significant 
(AZ k, 0.2 dex) only at metallicities 12+log(0/H) < 7.5 
and 12+log(0/H) > 8.5. Although the difference in the 
lowest-metallicity range is not a serious problem (because 
in this metallicity range the expected [Nn]A6584 flux is 
extremely weak and thus its measurement would be very 
challenging and probably inaccurate), it is important to 
pay attention to the difference in the high-metallicity 
range. Note that such "high-metallicity" range (where 
the discrepancy with Pettini & Pagel 2004 occurs) is 
not so metal rich — the metallicity 12+log(0/H) = 
8.5 corresponds to Z = 0.65Z Q , still in the sub-solar 
metallicity domain. 

In Figure 12, we compare the empirical calibrations 
of F([Oiii]A5007)/F([Nii]A6584 derived by us with the 
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one derived by Pettini & Pagel (2004). The difference be- 
tween the two calibrations is more serious than that seen 
in Figure 11 in the low metallicity range, 12+log(0/H) 
< 8. However Pettini & Pagel (2004) correctly mentioned 
that the flux ratio F([Oiii]A5007)/F([Nii]A6584 is of lit- 
tle use when F([Om]A5007)/F([Nii]A6584 > 100 because 
of the saturation of this diagnostic. The behavior of this 
diagnostic flux ratio in the low-metallicity range would be 
important to derive the upper limits on the metallicity 
from an upper limit of the [Nn]A6584 flux. 

4.3. Comparison with photoionization models 

To interpret the metallicity dependences of the emission- 
line flux ratios, we compare observational data with the 
predictions of photoionization models. In Figures 13 and 
14, we show the empirical metallicity dependences and 
the theoretical metallicity dependences of some metallic- 
ity diagnostics, where the latter are taken from Kewley 
& Dopita (2002) except for F([Nn]A6584)/ J F(Ha) that is 
taken from Kobulnicky & Kewley (2004). Since the ex- 
plicit analytic expression for the metallicity dependence 
of J F 1 ([Om]A5007)/F([Oii]A3727) is not given by Kewley 
& Dopita (2002), we derive the polynomial expression of 
the theoretical metallicity dependence by fitting the re- 
sults given in Table 2 of Kewley & Dopita (2002). The 
photoionization models presented by Kewley & Dopita 
(2002) and Kobulnicky & Kewley (2004) were calculated 
by the photoionization code MAPPINGS III (Sutherland 
& Dopita 1993) combined with the stellar population syn- 
thesis codes PEGASE (Fioc & Rocca-Volmerange 1997) 
and STARBURST99 (Leithcrcr et al. 1999), for the range 
7.6 < 12+log(0/H) < 9.4. They assume that stars and 
gas have the same metallicity, which is a reasonable as- 
sumption given that photoionization is due to hot, young 
stars, presumably recently formed from the same gas that 
they arc photoionizing. In their calculations, nitrogen is 
assumed to be a secondary nucleosynthesis element at 
12+log(0/H) > 8.3, and a primary nucleosynthesis ele- 
ment at lower metallicity. Effects of dust grains on the de- 
pletion of gas-phase heavy elements and on the radiative 
transfer are consistently taken into account. Their calcu- 
lations cover the range of ionization parameters —3.8 < 
log U < -2.0, or equivalently, 5 x 10 6 cm s _1 < q < 3 x 
10 8 cm s" 1 (where U = q/c). See Kewley & Dopita (2002) 
for details on the calculations. Note that they adopted 
12+log(O/H) = 8.93 (Anders & Grevesse 1989) and ex- 
pressed the metallicity in units of Z Q (12+log(O/H) = 
8.93). However, since we adopt a more recent value for the 
solar abundance, 12+log(O/H) = 8.69 (Allende Prieto et 
al. 2001), the metallicity notation is different when the Zq 
unit is used, which should be kept in mind to compare our 
results with their predictions. 

The most remarkable matter in the comparison be- 
tween the empirical and theoretical metallicity depen- 
dences of emission-line flux ratios is the significant dis- 
crepancy in the theoretically-expected i?23-sequence with 



respect to the observed trend. This is especially signifi- 
cant at low metallicity range 12+log(0/H) < 8. Shi ct 
al. (2006) also recently reported that a previous theoret- 
ical calibration of R23 (see McGaugh 1991; Kobulnicky 
et al. 1999) overpredicts the gas metallicity with respect 
to the metallicity measured through the gas temperature 
determined with [Om]A4363 line (AZ ~ 0.2 dex), espe- 
cially in the low metallicity range [i.e., 12+log(0/H) < 
8]. This discrepancy is not due to an improper compi- 
lation in our data, because it has been reported also in 
the earlier works that the empirical peak of R23 is seen 
around 12+(0/H) ~ 8.0, as mentioned already in §3. The 
discrepancy cannot be ascribed to problems to the model 
results of Kewley & Dopita (2002) either, because other 
theoretical works also predict higher peak metallicity of 
i?23 independently [12+log(0/H) 8.3; e.g., Kobulnicky 
et al. 1999]. One possible idea to reconcile this discrepancy 
is that the ionization parameter of the gas is higher than 
the parameter range which Kewley & Dopita (2002) cov- 
ers, especially in low-metallicity objects. If the ionization 
parameter correlates negatively with the gas metallicity 
and it reaches up to log U > — 2 at the lowest metallicities, 
photoionization models would predict larger values of R23 
in the lower-metallicity range with respect to constant- 
ly models. This idea appears to be consistent with the 
behaviors of the empirical sequences in the {/-sensitive 
flux ratios, F([Om]A5007)/F([Nn]A6584) (Figure 13) and 
F([Oiii]A5007)/ j F([Oii]A3727) (Figure 14). By focusing on 
these two [/-sensitive flux ratios, we can see that the ion- 
ization parameter increases by ~0.7 dex with decreasing 
oxygen abundance from 12+log(0/H) = 9.0 to 7.5, sup- 
porting the above interpretation. Although the absolute 
value of the required ionization parameter appear to be 
inconsistent between R23 and the latter two {/-sensitive 
flux ratios, the inferred absolute U values depends also 
on some model assumptions such as the spectral energy 
distribution (SED) of ionizing photons or the relative el- 
emental abundance ratios, which also change as a func- 
tion of metallicity. We thus conclude that the metallicity 
dependence of the ionization parameter (hereafter ll U-Z 
relation") causes the discrepancy between the empirical 
R23 distribution and the model predictions with a con- 
stant ionization parameter. 

Note that the F([Om]A5007)/F([Nll]A6584) and 
F([Oiii]A5007)/F([Oii]A3727) ratios are also sen- 
sitive to the hardness of the ionizing radiation, 
which is a strong function of the stellar metallic- 
ity. This effect can in principle also contribute to 
the dependence of F([Om]A5007)/F([Nll]A6584) and 
F([Om]A5007)/i ;l ([Oii]A3727) ratios on metallicity. 
However, the models by Kewley & Dopita (2002) plotted 
in Figures 13 and 14 already take into account the 
hardening of the stellar spectra as a function of metal- 
licity. Therefore, the discrepancy between constant-?/ 
models and the data indicates that the hardening of the 
ionizing spectra must be associated with a variation of 
U with metallicity. In particular, the dependence of the 
F([Om]A5007)/F([Oil]A3727) ratio on metallicity cannot 
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entirely be ascribed only to the hardening of the ionizing 
radiation, but also to a U-Z relation. 

The inferred U-Z relation is a very interesting result. 
Maier et al. (2006) also recently reported that the lower- 
metallicity galaxies tend to be characterized by a higher 
ionization parameter (see also Maier et al. 2004, who re- 
ported the correlation between the absolute B magnitude 
and the flux ratio of F([Om]A5007)/F([Oii]A3727) among 
galaxies in the local universe). Although a detailed theo- 
retical interpretation of this empirical relation goes be- 
yond the scope of this paper, in the following we dis- 
cuss two possible qualitative interpretations. One possi- 
ble origin of this effect may be associated with the mass- 
metallicity relation and with the mass-age relation in local 
galaxies. According to these relations, higher metallicity 
galaxies are associated with more massive and older sys- 
tems. Hn regions ionized by later stellar populations are 
expected to be characterized by lower ionization param- 
eters, due to the lower luminosity of the ionizing stars. 
Another possible explanation may be the (plausible) re- 
lation between gas metallicity and stellar metallicity, and 
in particular that lower metallicity gas is ionized by lower 
metallicity stars. For a given stellar mass, lower metallic- 
ity stars emit a harder and stronger radiation field, there- 
fore giving a higher ionization parameter. The latter ef- 
fect would naturally yield a U-Z relationship. The former 
are just qualitative interpretations. However, a thorough 
investigation of this phenomenon will requite detailed ob- 
servational studies of stellar population in star forming 
galaxies. 

The comparison of the empirical and the the- 
oretical sequences of the two [/-sensitive diagnos- 
tic flux ratios, F([Olll]A5007)/F([Nll]A6584) and 
F([Oiii]A5007)/F([Oii]A3727), also suggests the fact that 
the dispersion of the ionization parameter for a given 
metallicity should be relatively small. The typical RMS of 
the two flux ratios are ~0.5 (in logarithm) at 12+log(0/H) 
~ 7.5. This corresponds to an RMS of the ionization 
parameter of ^0.5 dex. This is the reason why the very 
[/-sensitive flux ratio, F([Oiii]A5007)/F([Oii]A3727), 
shows a clear metallicity dependence as seen in Figure 
13. The [/-metallicity relationship is also important to 
understand the behavior of the empirical metallicity de- 
pendence of the flux ratio F([Om]A5007)/F([Nll]A6584). 
This flux ratio is predicted to decrease with the oxygen 
abundance below 12+log(0/H) ~ 7.6 by photoionization 
models with a constant ionization parameter. Owing to 
the metallicity dependence of the ionization parameter, 
this flux ratio does not show the "turnover" seen in i? 23 
and thus it is very useful to investigate the gas metallicity 
of galaxies without the measurement of F([0iii]A4363). 
Another implication of these results is that one should 
not use constant-?/ photoionization models to derive the 
oxygen abundance from the observed flux ratios, not 
only from F([Om]A5007)/F([Nll]A6584 but also from any 
other metallicity diagnostics, which introduce systematic 
errors in the calibration. The empirical relations provided 
in this paper (Tables 5 and 6) are very useful to avoid 



such systematic errors to derive the gas metallicity by 
using only strong emission lines. 

As for the [/-insensitive diagnostic flux ratios, 
F([Nn]A6584)/F(Ha), F([Nii]A6584)/F([Oii]A3727) and 
F([Nll]A6584)/F([Sn]A6720), there arc no significant dis- 
crepancies between the empirical sequence and the the- 
oretical sequence (with a constant ionization parame- 
ter). This indirectly supports the above interpretation 
that the apparent discrepancy in i? 23 between the em- 
pirical sequence and the results of photoionization model 
is caused by the effect of the ionization parameter. 
Note that there is little or no metallicity dependence 
of the flux ratios of F([Nii]A6584)/F([Oii]A3727) and 
F([Nil]A6584)/ J F([Sll]A6720) in the low-mctallicity range, 
12+log(0/H) J5 8.0, in terms both of empirical and theo- 
retical dependences. Therefore these diagnostic flux ratios 
are useful only for the high metallicity galaxies. 

The photoionization models presented in Figures 13 
and 14 suggest an additional interpretation of the dis- 
crepancy in some diagnostics between the two samples 
discussed in §§4.1 (i.e., the discontinuity between sample 
A+B and sample C). Focusing on the metallicity range of 
12+log(0/H) ~ 8.3 where the two datasets of sample A+B 
and sample C overlap, we note that the trend of the dis- 
crepancy suggests that the galaxies in sample A+B have 
higher ionization parameter than the galaxies in sample C. 
This supports the interpretation that the discrepancy is 
at least partly caused by the selection effect, i.e., galaxies 
with higher ionization parameter are selectively picked up 
in sample A+B. Then, what causes this selection effect? 
This may be related with the fact that the [0m]A4363 
emission is extremely weak in higher metallicity galaxies. 
This means that we can measure the [Oiii]A4363 flux of 
galaxies with 12+log(0/H) ~ 8.3 (the highest metallic- 
ity in the galaxies in sample A+B) only when the [Oin] 
emission is very strong, which corresponds to a very high 
ionization parameter. 

4.4. Implications for studies of high-redshift galaxies 
and new diagnostics 

Although the i? 2 3 method is thought to be a good metal- 
licity diagnostic, various other diagnostics (some of which 
are investigated in this paper) have been proposed up to 
now. Indeed one of the main problems of the i? 23 method is 
that there are two solutions for a given i? 23 value and thus 
one cannot obtain a unique metallicity solution. Most of 
the newly proposed diagnostics use the [Nii]A6584 line to 
remove the degeneracy in i? 23 , because the secondary nu- 
cleosynthesis of nitrogen makes this line emission very sen- 
sitive to the gas metallicity. However, there are two non- 
negligible problems with the use of the [Nn]A6584 line. 
First, especially for low-metallicity systems, the contri- 
bution of the primary nucleosynthesis and the secondary 
nucleosynthesis in the nitrogen abundance is not well un- 
derstood, which leads to an uncertainty in the relative ni- 
trogen abundance as a function of the metallicity. Second, 
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the [Nn]A6584 emission is in the red part of the rest-frame 
optical spectrum of galaxies, which prevents its applica- 
tion to the observational investigations of high-z systems. 
For example, the optical detectors with a sensitivity up to 
A ~ l^m can detect the [Nn]A6584 emission of galaxies 
only at z < 0.52, and the F-band atmospheric window 
limits the highest redshift to z ~ 2.7 for ground-based fa- 
cilities. Although one of the undoubtfully interesting tar- 
gets for the JWST is the population related to the cos- 
mic reionization, the sensitivity of NIRSpec (Posselt et 
al. 2004) boarded on JWST can examine the [Nn]A6584 
emission of the objects at z SS 6.6, where the cosmic reion- 
ization has already nearly ended (e.g., Kashikawa et al. 
2006; Fan et al. 2006). Another problem associated with 
the [Nn]A6584 line is that it becomes very weak and diffi- 
cult to measure at low metallicities: [Nn]A6584/Ho: < 0.1 
at 12+log(0/H) < 8.5. 

Our results on the empirical metallicity depen- 
dences suggest that one does not need [Nn]A6584 any 
more to distinguish the upper- and lower-branches 
of the i?23 sequence. This is because the flux ra- 
tio of F([Oin]A5007)/F([Oii]A3727) is also a good 
metallicity diagnostics, thanks to the small disper- 
sion of the ionization parameter at a given metallic- 
ity. The empirical R23 sequence peaks at 12+log(0/H) 
~ 8.0, where the empirically determined flux ratio of 
F([Om]A5007)/F([Oli]A3727) is -2. Therefore one can 
recognize whether the observed R23 belongs to the upper- 
branch of the i?23 sequence or not, depending on whether 
F([Oiii]A5007)/F([Oii]A3727) < 2 or not. Note that this 
result is consistent with an earlier remark by Maier et al. 
(2004) that the flux ratio of F([Oiii]A5007)/F([Oii]A3727) 
can be used to distinguish the upper- and lower-branches 
of the i?23 sequence. Our work gives the physical expla- 
nation for this idea (the U-Z relation) and a criterion to 
distinguish the degeneracy [F([Om]A5007)/F([Oii] A3727) 
< 2] on the remark by Maier et al. (2004). 

The above result is due to the fact that the ion- 
ization parameter has a strong metallicity dependence, 
and it thus implies that the ionization parameter itself 
is a sort of metallicity diagnostic. Motivated by this, 
we examine the metallicity dependence of the flux ratio 
F([Neiii]A3869)/F([On]A3727), in Figure 15. The reasons 
for focusing on this flux ratio are: (a) the two emission 
lines have different ionization degrees, their ratio should 
have a strong dependence on the ionization parameter and 
therefore is a possible good metallicity diagnostics; (b) 
their wavelength separation is very small and thus their 
flux ratio is not significantly affected by dust reddening; 
and (c) the two lines are located at a blue part in the rest- 
frame optical spectrum and thus their flux ratio could be a 
powerful diagnostic even for high- z galaxies. As expected, 
this flux ratio shows a clear metallicity dependence, which 
is apparently seen in Figure 14. In Tables 7 and 8, the 
mean and the RMS of this flux ratio for within each bins 
of oxygen abundance are given, just similar to Tables 3 
and 4 (§3). To obtain the analytic expression of this re- 
lation, we fit the observed sequence with a second-order 



polynomial function. The coefficients of the fit are given 
in Tables 9 and 10. 

This flux ratio can be measured for galaxies up to 
z <~ 1.6 with optical instruments, up to z <~ 5.2 with 
near-infrared instruments on the ground-based facilities, 
and up to z ~ 12 with JWST/NIRSpec, therefore this 
flux ratio is a promising tool for metallicity studies at 
high redshift. In particular, it is useful for low metallicity 
galaxies, for which the intensity of [Nein]A3869 becomes 
comparable to [On]A3727 and therefore easier to detect 
[F([Nem]A3869)/F([Oii]A3727) > 0.2 at 12+logO/H <8]. 
Detailed theoretical calibrations on this flux ratio are re- 
quired, taking the metallicity dependence of the ionization 
parameter into account, which go beyond the scope of this 
paper. 

One possible caveat for the use of the diagnos- 
tic flux ratios of F([Nem]A3869)/F([Oii]A3727) [and 
F([Oiii]A5007)/F([Oii]A3727), too] may be the ef- 
fect of AGNs. Since AGNs also tend to show 
higher ratios of F([Nein]A3869)/F([Oii]A3727) and 
F([Offl]A5007)/F([Oll]A3727), galaxies harboring an 
AGN may be misidentified as low-metallicity galax- 
ies. However, we can identify AGNs through the de- 
tection of HeiiA4686 and/or [Nev]A3426. Nagao et 
al. (2001) reported that typical type-2 AGNs show 
F([Nev]A3426)/F([On]A3727) -0.4, and typical type-1 
AGNs show even higher ratio (<; 1). Lamareille et al. 
(2004) also reported that AGNs and star-forming galaxies 
can be distinguished by using diagnostic diagrams using 
only the blue part of the spectrum, i.e., O32 versus F23 
and F([Oni]A5007)/F(H/?) versus F([Oii]A3727)/F(H/3) 
(see also Rola et al. 1997). These suggest that we can 
easily distinguish AGNs from low-metallicity galaxies by 
using only diagnostics available in the blue part of the 
spectrum, even with moderate quality spectroscopic data. 
Another caveat for the use of some diagnostic flux ratios 
calibrated in this paper especially for high-z galaxies is 
that several of the empirical relations rely on the U-Z re- 
lation. It is not obvious that the U-Z relation found in the 
local galaxies also holds for high-z galaxies. If the U-Z re- 
lation is a consequence of the relation between gas and 
stellar metallicity, as discussed in the previous section, 
then the relation is not expected to evolve and should re- 
main valid at any redshift. Instead, if the U-Z relation is a 
consequence of the mass-metallicity relation which evolves 
with redshift (Savaglio et al. 2005; Erb et al. 2006; see also 
Maier et al. 2004), then also the U-Z relation may evolve 
with redshift and may require a re-calibration of our em- 
pirical relations at high redshift. The latter case would 
be a serious problem for several studies at high redshift. 
Indeed, most of the gas metallicity diagnostics discussed 
in this paper, including the ones most widely used (e.g. 
F23), are significantly affected by the dependence on the 
ionization parameter. 

Another difficulty to measure the gas metallicities of 
high-z galaxies is the faintness of targets, which some- 
times prevents from measuring accurate emission-line 
fluxes. The use of low- resolution grating to improve the 
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signal-to-noise ratio may yield to a blending of the Ha 
and [Nil] emission lines, which results in poor deter- 
minations of the gas metallicity. Therefore It may be 
useful to investigate metallicity diagnostics which use 
the sum of F(Ha) and F([Nn]). In particular, we have 
examined the metallicity dependence of the flux ra- 
tio F(Ha+[Nn]AA6548,6584)/F([Sn]A6720) in Figure 16. 
This group of lines can be measured even in low-resolution 
spectra and even in spectra covering a relatively narrow 
wavelength range, and therefore may be particularly use- 
ful in high-z studies. There is a clear dependence of this 
flux ratio on the oxygen abundance, seen as a U-shaped 
distribution with a minimum at 12+log(0/H) <~ 8.7 (i.e., 
■^gas ~ Zq). The mean and the RMS of this flux ratio for 
each bin of oxygen abundance are given in Tables 7 and 
8, and the coefficients of the fit are given in Tables 9 and 
10. The observed distribution of this flux ratio is naturally 
expected, since the behavior of the nitrogen emission as 
a secondary element should dominate at the super-solar 
metallicity range, while F([Nn]) and F([Sn]) should be- 
come weak with respect to F(Ha) at the low-metallicities 
due to the decrease of the corresponding ions. Although 
this diagnostic like i?23 has two solutions when the ratio is 
below 10, this ratio seems useful for low-metallicity galax- 
ies where it is larger than 10, in which case it is possible 
to state that the object belongs to the lower branch of 
the U-shaped distribution. This diagnostic is also useful 
when the [Sn] emission is not detected (this is frequently 
the case when high-z faint galaxies are concerned). In this 
case, we can calculate a lower limit for this flux ratio, and 
we can derive accordingly an upper limit to the gas metal- 
licity if the lower limit is larger than 10. Note that this 
diagnostic is essentially independent of dust reddening. 

Finally, We have also investigated the flux ra- 
tios of F([Om]A5007)/F(H/?) and F([On]A3727)/F(H/3) 
(Figures 17 and 18). The empirical calibrations for these 
two flux ratios may be useful when either [On] or [Om] 
are not available, because out of the wavelength range or 
on a strong OH airglow emission line, or in a region of 
bad atmospheric transmission. The means and the RMSs 
of these two flux ratios for each bin of oxygen abundance 
are given in Tables 7 and 8, and the coefficients of the 
fit are given in Tables 9 and 10. As expected, both of 
the two flux ratios again show fl-shaped distributions just 
similar to the R 23 parameter. The F([Om]A5007)/F(H/3) 
ratio is useful especially for high metallicity galaxies, be- 
cause the targets should belong to the upper-branch of the 
distribution when F([Om]A5007)/F(H/3) < 1 (although 
this might be wrong when extremely metal-poor galaxies 
[12+log(0/H) < 7.0] are concerned). Note that this flux 
ratio is very sensitive to the oxygen abundance and the 
dispersion of the data is small at F([Om]A5007)/F(H/3) 
< 1. 

Figure 19 summarizes the use of some of the metallic- 
ity diagnostics discussed in this paper as a function of red- 
shift and for various observing facilities, and in particular 
optical spectrometers, ground-based near-IR spectrome- 
ters and NIRSpec on board of JWST. In principle (i.e., 



sensitivity permitting), MIRI on board of JWST will be 
able to observe the same diagnostics at even higher red- 
shifts. Note that the ratio F([Neill]A3869)/F([On]A3727) 
extends the diagnostic capability of any observing facility 
to significantly higher redshift. 

5. Summary 

We have combined two large spectroscopic datasets to 
derive empirical calibrations for gas metallicity diagnos- 
tics involving strong emission lines. The two datasets con- 
sist of about 50000 spectra from the SDSS DR4, which 
probe metallicities 12+log(0/H)>8.3 (sample C), and of 
328 spectra of low metallicity galaxies with a measure- 
ment of the [Oiii]A4363 line (sample A+B), which probe 
metallicities 12+log(0/H)<8.4. Together, these two sam- 
ples provide the largest dataset of galaxies with known 
metallicity currently available, and spanning more than 2 
dex in metallicity. 

We have provided empirical calibrations both for 
metallicity diagnostics already proposed in the past and 
for new metallicity indicators proposed in this paper. We 
have given an analytical description for the metallicity de- 
pendence of the following diagnostics and line ratios: F23, 
F( [Nil] A6584) /F (Ha) , F( [Om] A5007) /F ( [Nil] A6584) , 
F([Nii]A6584)/F([Oii]A3727), F([Nn]A6584)/ 
F([Sn]A6720), F([Oiii]A5007)/F([Oii]A3727), and 
F([Nelll]A3869)/F([Oll]A3727). The calibrations are per- 
formed within the metallicity range 7.0 < 12+log(0/H) 
< 9.2. All of the investigated flux ratios show strong 
dependences on metallicity, at least in some metallicity 
ranges. We have shown that the monotonic metallicity 
dependence of the ratio F([Oiii]A5007)/F([Oii]A3727) 
can be used to break the degeneracy of the F23 
parameter when F([Nn]A6584)/F(Ha) is not avail- 
able. The F([Oiii]A5007)/F([Oii]A3727) ratio is 
particularly useful at high redshift, where Ha and 
[Nn]A6584 are shifted outside the observed band. 
Another promising metallicity tracer at high-z 
is the ratio F([Neiii]A3869)/F([Oii]A3727), which 
is found to anti-correlate with metallicity. The 
F([Neiii]A3869)/F([On]A3727), ratio is particularly 
useful at high redshift, where most of the other diagnostic 
lines are shifted outside the observed band. 

We have also investigated the observed relationships 
through a comparison with photoionization models. Some 
of the diagnostics investigated in this paper are strongly 
dependent on the ionization parameter U. The observed 
trends of these diagnostics highlight a clear, inverse rela- 
tionship between ionization parameter and metallicity in 
galaxies. Such a strong U-Z relationship is also required 
to explain the trend observed for the F23 parameter. The 
U-Z relationship is relatively tight and, indeed, we have 
found that at any given metallicity the ionization param- 
eter has a small dispersion (~0.5 dex). The strong re- 
lationship between ionization parameter and metallicity 
in galaxies should warn about the use of simple models, 
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which assume constant ionization parameter, to infer gas 
metallicities from line ratios. 
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Fig. 2. Emission-line flux ratios of [Sii]A6717/[Sii]A6731 
for the compiled low- metallicity galaxies, as a function of 
the oxygen abundance derived by us (see §§2.1.2). The 
upper horizontal dotted line denotes the theoretical low- 
density limit of this flux ratio and the lower dotted line 
denotes the high-density limit. 



O 

cn 
o 




7.5 



7 

7 7.5 8 8.5 

12 + loq(0/H),, , 

=> v ' / hterature 

Fig. 3. Oxygen abundances of the compiled low- 
metallicity galaxies re-calculated by us (see §§2.1.2) are 
plotted as a function of the oxygen abundances given in 
the original references. Dotted line is not the best-fit line 
but a reference line for the case when the two quantities 
are the same. 
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Fig. 4. Frequency distribution of the redshift of the SDSS 
DR4 galaxies after our sample selection (sample C) de- 
scribed in §§2.2. Galaxies at z < 0.028 are not included 
(see text). 
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Fig. 6. Emission-line flux ratios of R23 (= [F( [On] A3727) Fig. 7. Same as Figure 6 but for the emission- 

+ 1.327xF([Oin]A5007)]/F(H ; 3)), F([Nii]A6584)/F(Ha), line flux ratios of F([Nll]A6584)/F([Oll]A3727), 

and F([Oiii]A5007)/F([Nii]A6584) for galaxies in sample F([Nll]A6584)/F([Sll]A6720), and 

A (red triangles), in sample B (green circles) and in sample F([Oni] A5007)/i 7 '([Oii]A3727). 

C (black dots), as a function of the oxygen abundance. 

The compiled low-metallicity galaxies with an error of the 

oxygen abundance larger than 0.05 dex are not plotted. 

Dotted lines denote the solar metallicity [12+log(0/H) = 

8.69]. 
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Fig. 9. Same as Figure 8 but for the emission- 
line flux ratios of F([Nii]A6584)/F([Oii]A3727), 
F([Nii]A6584)/F([Sii]A6720), and 
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Fig. 10. Comparison of our results with the previous 
empirical metallicity calibrations for the R23 parameter. 
Solid red line denotes our calibration. Blue, green, and 
magenta lines denote the calibration given by Tremonti et 
al. (2004), Edmunds & Pagel (1984), and Zaritsky et al. 
(1994), respectively. Symbols and errorbars are the same 
as those in Figure 8. Vertical dotted line denotes the solar 
metallicity [12+log(0/H) = 8.69]. 
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Fig. 11. Comparison of our results with the previous em- 
pirical metallicity calibrations for F([Nn]A6584)/F(Ha). 
Solid red line denotes our calibration. Blue and magenta 
lines denote the calibration given by Pettini & Pagel 
(2004) and Denicolo et al. (2002), respectively. Symbols 
and errorbars are the same as those in Figure 8. Vertical 
dotted line denotes the solar metallicity [12+log(0/H) = 
8.69]. 
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Fig. 12. Comparison of our results with the 
previous empirical metallicity calibrations for 
F([Oiii]A5007)/F([Nii]A6584). Solid red line denotes 
our calibration, and blue line denotes the calibration 
given by Pettini & Pagel (2004). Symbols and errorbars 
are the same as those in Figure 8. Vertical dotted line 
denotes the solar metallicity [12+log(0/H) = 8.69]. 
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Fig. 13. The averaged flux ratios and the best-fit poly- Fig. 14. Same as Figure 13 but for the emission- 

nomial functions of the metallicity dependence of i? 2 3, line flux ratios of F([Nii]A6584)/F([Oii]A3727), 

J F([Nii]A6584)/F(Ha) and F([Oiii]A5007)/F([Nii]A6584) F([Nii]A6584)/F([Sii]A6720), and 

(dashed line) are compared with the predictions of pho- F([Om]A5007)/i ;l ([Oii]A3727). 

toionization models (solid lines: Kewley & Dopita 2002; 

Kobulnicky & Kewley 2004). The lines with a digit 1, 2, 

3, 4, 5, 6, and 7 denote the model predictions with the 

ionization parameter of log U — -3.8, -3.5, -3.2, -2.9, - 

2.6, -2.3, and -2.0, respectively. Dotted line denotes the 

solar metallicity [12+log(0/H) = 8.69]. 
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Fig. 15. (Upper) Emission-line flux ratios of Fig. 16. Same as Figure 15 but for the emission-line flux 

F([Nem]A3869)/ F([On]A3727) of the galaxies in sample ratio of F(Ha+[Nii]AA6548,6584)/F([Sii]A6720). 

A (red triangles), those in sample B (green circles), and 

those in sample C (black dots), as a function of the 

oxygen abundance. As for sample C, only the objects 

with S/N([Nem]) > 10 and S/N([0n]) > 10 (cataloged 

values) are plotted. (Lower) Same as the upper panel 

but the mean and the RMS values are shown for each bin 

of the oxygen abundance, instead of the individual data. 

Filled stars and filled squares denote the mean flux ratios 

for galaxies in sample A+B and those for galaxies in 

sample C, respectively. The errorbar denotes the RMS for 

each metallicity bin. The dashed line denotes the best-fit 

polynomial (second-order) function. 
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Fig. 17. Same as Figure 15 but for the emission-line flux Fig. 18. Same as Figure 15 but for the emission-line flux 
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Fig. 19. Schematic view of the availability of various metallicity diagnostics for each redshift. The black solid curves 
indicate the effect of redshift for some of the diagnostic lines discussed in this paper. The colored boxes indicate 
the wavelength coverages of optical spectrometers (blue), of ground-based near-IR spectrometers (magenta), and of 
NIRSpcc/MIRI on board of JWST (red). The marks on the right of the diagram indicate the maximum redshift at 
which some of the metallicity diagnostics can be used with the various facilities. 
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Table 1. Re-calculated properties of the compiled low- metallicity galaxies 



Object R 23 ™ H (S+) a i(0 2+ ) b 12 + log(0/H) Ref. c 



HS 0029+1748 


8.624 


± 


0.168 


70+ 47 

' u -42 


1.286 + 0.016 


8. 


.046 


± 


0. 


016 


I04b 


HS 0111+2115 


9.142 


± 


0.179 


< 22 


1 inc+0.063 

i.ius_ 066 


8. 


290 


± 


0. 


076 


I04b 


HS 0122+0743 


6.669 


± 


0.127 




1.777 ± 0.025 


7. 


.597 


± 


0. 


014 


I04b 


HS 0128+2832 


10.272 


■ ± 


0.197 


104tf 7 


1.256 + 0.011 


8. 


.147 


± 


0. 


013 


I04b 


HS 0134+3415 


10.485 


■ ± 


0.204 


200+™ 


1.639 + 0.019 


7. 


.858 


± 


0. 


014 


I04b 


HS 0735+3512 


9.640 


± 


0.163 


7C+39 
'°-36 


1.205 ± 0.015 


8. 


.193 


± 


0. 


014 


I04b 


HS 0811+4913 


9.861 


± 


0.184 


44+^1 


1.449 + 0.016 


7. 


.968 


± 


0. 


013 


I04b 


HS 0837+4717 


8.061 


± 


0.153 


373±} I 


1.954 ± 0.026 


7. 


587 


± 


0. 


013 


I04b 


HS 0924+3821 


8.617 


± 


0.152 


16±*| 


1.256 ± 0.021 


8. 


.089 


± 


0. 


018 


I04b 


HS 1028+3843 


10.464 


± 


0.206 


492+^5 


1.582 + 0.017 


7. 


.891 


± 


0. 


014 


I04b 


HS 1213+3636A 


7.353 


± 


0.112 


or +35 
Jd -32 


1 074+ ' 026 

1 - u '^-0.027 


8. 


263 


± 


0. 


033 


I04b 


HS 1214+3801 


9.095 


± 


0.162 


20+ 37 


1.339 + 0.013 


8. 


026 


± 


0. 


012 


I04b 


HS 1311+3628 


8.570 


± 


0.146 


951^ 


1.141+0.013 


8. 


199 


± 


0. 


015 


I04b 


HS 2236+1344 


6.975 


± 


0.135 


86^ 


2.123 + 0.032 


7. 


.464 


± 


0. 


013 


I04b 


HS 2359+1659 


9.608 


± 


0.125 


< 10 


1 18Q+ ' 016 


8. 


179 


± 


0. 


017 


I04b 


IC 0010 1 


6.757 


± 


0.400 


(10) d 


1 01 c,+ 077 
i.uio_ 0g5 


8. 


301 


± 


0. 


096 


L03b 


IC 0010 2 


6.607 


± 


0.369 


(10) d 


QS4+ ' 077 


8. 


.335 


± 


0. 


097 


L03b 


IC 0010 3 


7.658 


± 


0.624 


(10) d 


1 236+ ' 155 


8, 


083 


± 


0. 


128 


L03b 


IC 1613 


7.641 


± 


0.822 


< 167 


1 796+ ' 161 
r. '»o_ 15g 


7. 


.643 


± 


0. 


085 


L03a 


IC 5152 


5.852 


± 


0.509 


91+275 
zi -21 


1 240+ ' 166 

4- zw -0.183 


7. 


955 


± 


0. 


135 


L03a 


KISSB 0023 


5.251 


± 


0.238 


< 172 


1 OQQ+0.094 


7. 


.570 


± 


0. 


030 


M04 


KISSB 0061 


7.348 


± 


0.239 


< 10 


1 565+ 037 

J--JOJ-0 .038 


7. 


.772 


± 


0. 


024 


L04 


KISSB 0086 


8.177 


± 


0.242 


< 10 


1.241 +0.028 


8. 


.088 


± 


0. 


025 


L04 


KISSB 0171 


8.276 


± 


0.247 


< 35 


1.189 + 0.020 


8. 


121 


± 


0. 


022 


L04 


KISSB 0175 


9.346 


± 


0.295 


166ll» 


1.344 ± 0.024 


8. 


.041 


± 


0. 


021 


L04 


KISSR 0049 


7.786 


± 


0.398 




1 Qne+0.083 
1 - OUO -0.086 


8, 


.028 


± 


0. 


059 


M04 


KISSR 0073 


7.424 


± 


0.382 




1 Sfifi+O- 056 
l.JDD_ 057 


7. 


.943 


± 


0. 


042 


L04 


KISSR 0085 


5.316 


± 


0.222 


oqq+727 


1 77K+0.104 
' d -0.103 


7. 


.532 


± 


0. 


043 


M04 


KISSR 0087 


7.827 


± 


0.247 


71 +iii 

' L -71 


0.996 ± 0.023 


8. 


.386 


± 


0. 


030 


L04 


KISSR 0116 


7.826 


± 


0.241 


< 90 


1 9ne;+ ' 021 


8. 


.107 


± 


0. 


022 


L04 


KISSR 0286 


7.710 


± 


0.235 


29ti° 3 


1.103 + 0.022 


8. 


.216 


± 


0. 


028 


L04 


KISSR 0310 


9.501 


± 


0.551 


9R +198 
zo _26 


1.528 + 0.050 


7. 


892 


± 


0. 


039 


L04 


KISSR 0311 


8.935 


± 


0.468 


q 9 +192 

oz -32 


1 007+0. 045 

l.OO (_Q Q 46 


7. 


.994 


± 


0. 


036 


L04 


KISSR 0396 


8.023 


± 


0.238 


27+J? 


1 4n6 +0 046 

1 ' 4uu -0 .047 


7. 


938 


± 


0. 


029 


M04 


KISSR 0666 


8.630 


± 


0.397 


< 10 


9 ! co+0.092 


7. 


.528 


± 


0. 


038 


M04 


KISSR 0675 


8.773 


± 


0.583 


1138+ 2735 


1 507+ ' 112 


7. 


890 


± 


0. 


071 


M04 


KISSR 0814 


8.620 


± 


0.276 


79+^f 


1.373 ± 0.026 


7. 


.983 


± 


0. 


022 


L04 


KISSR 1013 


7.271 


± 


0.295 


474+Jf 9 


1 781+0-105 
t- '81_0.104 


7. 


690 


± 


0. 


038 


M04 


KISSR 1194 


8.836 


± 


0.356 


ci +240 


1 /IKC + "-!' 41 

1.458I 042 


7. 


930 


± 


0. 


029 


M04 


KISSR 1490 


6.676 


± 


0.308 


< 381 




7. 


.568 


± 


0. 


040 


M04 


KISSR 1778 


6.628 


± 


0.325 


f;q+383 


i.ouo_ og4 


7. 


.955 


± 


0. 


064 


M04 


KISSR 1845 


9.440 


± 


0.381 


50± 240 


1.325 ± 0.033 


8. 


069 


± 


0. 


028 


M04 


Mrk 0005 


7.192 


± 


0.106 




1 220+ 051 

1 - zzu -0.052 


8. 


.058 


± 


0. 


041 


198 


Mrk 0022 


8.726 


± 


0.101 


7, +55 
'-■--49 


1 Q 4 g+0 020 
1 - J * a -0 .021 


8. 


002 


± 


0. 


.015 


194 


Mrk 0035 


7.892 


± 


0.121 


189+.H 


1 021+ 012 

± - uzl -0.013 


8. 


.368 


± 


0. 


015 


I04b 


Mrk 0036 


7.708 


± 


0.092 


< 102 


1.524 + 0.037 


7. 


.816 


± 


0. 


021 


198 


Mrk 0067 


9.309 


± 


0.176 


< 10 


1.320 ± 0.024 


8. 


059 


± 


0. 


019 


I04b 


Mrk 0116 


2.937 


± 


0.017 


86^ 6 2 


1.927 ± 0.038 


7. 


178 


± 


0. 


014 


197 


Mrk 0116 1 


2.935 


± 


0.043 


< 1369 


9 1 oq+0.064 


7. 


.084 


± 


0. 


023 


P92 


Mrk Olio 2 


3.009 


± 


0.071 


cq+722 
68 -68 


1 Q7Q+ ' 092 

r.y '»_o.o9o 


7. 


.217 


± 


0. 


027 


F92 


Mrk 0162 


8.180 


± 


0.083 


< 10 


1 104+0.043 

r- 1 »1_0.044 


8. 


.116 


± 


0. 


034 


198 


Mrk 0178 


8.517 


± 


0.245 


1 99 +234 
J.^-122 


1.588 + 0.104 


7. 


.816 


± 


0. 


057 


GOO 


Mrk 0193 


8.906 


± 


0.101 


172+H 


1.639 + 0.019 


7. 


.795 


± 


0. 


011 


194 


Mrk 0209 


8.075 


± 


0.018 




1.630 ± 0.007 


7. 


.755 


± 


0. 


.004 


197 


Mrk 0450 1 


8.514 


± 


0.144 


132t» 


1.173 + 0.013 


8. 


154 


± 


0. 


014 


I04b 


Mrk 0450 2 


8.661 


± 


0.166 


< 21 


1 9c;, +0.028 


8. 


094 


± 


0. 


.024 


I04b 


Mrk 0475 


8.392 


± 


0.111 


< 45 


1.411 + 0.028 


7. 


933 


± 


0. 


019 


194 
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Object 


R23 


n H (S+) a 


t(0 2 +) b 


12 + log(0/H) 


Ref. c 


Mrk 0487 


8.413 


± 


0.157 


— Do 


1.266+2,^ 

— U.058 


8.076 


± 


0.043 


197 


Mrk 0600 


8.578 


± 


0.101 


58+« 

— 38 


1.579 + 0.020 


7.824 


± 


0.012 


198 


Mrk 0724 


8.618 


± 


0.149 


19 + ?I 

^—19 


1.296+°°" 

^-0.015 


8.045 


± 


0.013 


I04b 


Mrk 0750 


8.357 


± 


0.079 


< 10 


1.205+°°^ 


8.128 


± 


0.021 


198 


Mrk 0930 


7.905 


± 


0.084 


56+12 

— 36 


1 ooc+0.037 
ou — 0.038 


8.084 


± 


0.029 


198 


Mrk 1063 


6.183 


± 


0.106 


96+« 

— 39 


1 027+22«2 

±.vj^ 1 _q.062 


8.260 


± 


0.069 


I04b 


Mrk 1089 


5.554 


± 


0.057 


92+« 

— 39 


1.108+°°?! 

— U.074 


8.090 


± 


0.070 


198 


Mrk 1236 


9.550 


± 


0.170 




1.225 + 0.012 


8.157 


± 


0.013 


I04b 


Mrk 1271 


9.680 


± 


0.074 


_ 51 


1.411 ± 0.018 


7.996 


± 


0.012 


198 


Mrk 1315 


9.104 


± 


0.164 


11+?? 

—11 


1.103 + 0.009 


8.270 


± 


0.013 


I04b 


Mrk 1328 


6.981 


± 


0.165 


25+?? 


0.937+2??2. 

w.utj 1 _q 219 


8.457 


± 


0.106 


V03 


Mrk 1329 


8.539 


± 


0.150 


18+?i 


1.080 + 0.009 


8.278 


± 


0.013 


I04b 


Mrk 1409 


8.754 


± 


0.138 


599+i?i 


1.362+2-222 

^^ — 0.067 


8.025 


± 


0.040 


197 


Mrk 1416 


8.098 


± 


0.065 


< 10 


1.514 + 0.031 


7.854 


± 


0.016 


197 


Mrk 1434 


7.640 


± 


0.049 


< 10 


1.551 ±0.015 


7.786 


± 


0.009 


197 


Mrk 1450 


7.669 


± 


0.052 


< 43 


1.330 + 0.016 


7.963 


± 


0.012 


194 


Mrk 1486 


8.140 


± 


0.059 


27+12 

— 27 


1.468 + 0.022 


7.884 


± 


0.013 


197 


NGC 2363 A 


9.358 


± 


0.020 


85+?* 

— 53 


1.584 + 0.006 


7.843 


± 


0.004 


197 


NGC 2363 B 


7.286 


± 


0.074 


14+™ 

— 14 


1.496+°°34 

— 0.034 


7.818 


± 


0.019 


197 


NGC 3109 


6.221 


± 


0.437 


(10) d 


1.463+2?2? 

j-.-iwfj — 0.277 


7.792 


± 


0.138 


L03b 


NGC 4214 A6 


6.898 


± 


0.162 


33 +93 


w ^— 0.050 


8.280 


± 


0.062 


K96 


NGC 4214 C6 


7.913 


± 


0.176 


74+1 

' — 74 


0.983+2^ 

— 0.02b 


8.426 


± 


0.029 


K96 


NGC 4861 


8.801 


± 


0.020 


74+^ 

1 —26 


1.363 + 0.006 


7.987 


± 


0.005 


197 


PGC 18096 


10.671 


± 


0.124 


195+21 

it -"-' — 57 


1.339 ± 0.021 


8.086 


± 


0.017 


GOO 


PGC 27864 1 


8.410 


± 


0.151 


116+1" 

—43 


1.648 ± 0.018 


7.771 


± 


0.012 


I04b 


PGC 27864 2 


7.616 


± 


0.135 


< 10 


1.657 + 0.027 


7.738 


± 


0.014 


I04b 


PGC 37727 


8.082 


± 


0.129 


57+22. 

— 36 


1.256+2'2?2 

^^ — 0.030 


8.079 


± 


0.022 


I04b 


PGC 39188 


7.325 


± 


0.091 


1088 + !n7 

uuv -107 


1.014+2-222 

w — 0.027 


8.375 


± 


0.030 


V03 


PGC 39402 


7.368 


± 


0.146 


261+Tnn 

— 100 


2.002 + 0.028 


7.518 


± 


0.014 


104 


PGC 39845 


5.990 


± 


0.181 


140+222 

w — 140 


1.689 + 0.146 


7.653 


± 


0.051 


V03 


PGC 40521 


5.549 


± 


0.350 


57+2,f 

— 57 


1.330+2?,! 

— 0.241 


7.864 


± 


0.145 


V03 


PGC 40582 1 


5.814 


± 


0.113 


108+™ 

ww — 66 


1.880 + 0.031 


7.480 


± 


0.015 


I04a 


PGC 40582 2 


5.033 


± 


0.103 


< 10 


1.829+22 5 J 

u "-0.056 


7.456 


± 


0.025 


I04a 


PGC 40582 3 


4.901 


± 


0.187 


< 62 


1.908+2^ 

— 0.32 l 


7.410 


± 


0.125 


I04a 


PGC 40582 4 


5.670 


± 


0.131 


(10) d 


1.918+2-212 

w — 0.077 


7.458 


± 


0.032 


I04a 


PGC 40604 


7.671 


± 


0.266 


< 127 


1.266 + S"i?2 

— 0. 175 


8.058 


± 


0.120 


V03 


PGC 40604 a 


6.668 


± 


0.279 


155+?S 

— 1 O 


1.275+Q- 1 " 

— 0.1 < 


7.972 


± 


0.118 


V03 


PGC 41360 


7.255 


± 


0.187 


39+ 107 

— 39 


1.542 + 0.047 


7.775 


± 


0.029 


V03 


PGC 42160 


5.418 


± 


0.238 


< 123 


1.486+°-^ 


7.744 


± 


0.104 


V03 


PGC 49050 


7.773 


± 


0.538 


< 157 


i.io3±g:i2l 


8.250 


± 


0.194 


L03a 


SBS 0335-052 


4.343 


± 


0.041 


275li?2 


2.040 + 0.036 


7.280 


± 


0.014 


198 


SBS 0335-052 E3 


4.588 


± 


0.090 


(10) d 


2.027 + 0.029 


7.306 


± 


0.014 


P06 


SBS 0335-052 E4-5 


4.443 


± 


0.093 


(10) d 


2.128l°-°S 


7.248 


± 


0.017 


P06 


SBS 0335-052 E7 


3.457 


± 


0.074 


(10) d 


1 074+O.O63 
0.062 


7.209 


± 


0.027 


P06 


SBS 0335-052 E-NW 


3.446 


± 


0.097 


(10) d 


2.001+2 loo 

— 0.122 


7.194 


± 


0.049 


P06 


SBS 0335-052 E-SE 


4.087 


± 


0.085 


(10) d 


1.979 + 0.053 


7.275 


± 


0.023 


P06 


SBS 0335-052 W 


2.550 


± 


0.072 


(10) d 


1 074+026O 

— 0.249 


7.133 


± 


0.073 


P06 


SBS 0749+568 


8.144 


± 


0.227 


< 10 


1.528 + 0.081 


7.843 


± 


0.045 


197 


SBS 0749+582 


11.727 


± 


0.272 


' —117 


1.334 + 0.034 


8.135 


± 


0.028 


197 


SBS 0907+543 


10.014 


± 


0.241 


118+?TI 


1.444 + 0.043 


7.974 


± 


0.031 


197 


SBS 0926+606 


8.117 


± 


0.059 


188+2 


1.434 + 0.023 


7.911 


± 


0.014 


197 


SBS 0940+544 


5.853 


± 


0.081 




2.016 + 0.038 


7.430 


± 


0.015 


197 


SBS 0943+561 


9.024 


± 


0.424 


971 +630 
z ' i -271 


L - '°°-0.129 


7.749 


± 


0.060 


197 


SBS 0948+532 


8.841 


± 


0.144 


7Q+85 
'°-73 


1.339 + 0.027 


8.014 


± 


0.021 


194 


SBS 1054+365 


8.959 


± 


0.090 


< 27 


1.383 + 0.019 


7.978 


± 


0.014 


197 


SBS 1116+583B 


7.014 


± 


0.249 


593+2JI 


1.670 ± 0.089 


7.673 


± 


0.049 


197 


SBS 1128+573 


8.570 


± 


0.212 


91 1 +339 

zii -211 


1.689 + 0.062 


7.751 


± 


0.033 


197 
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Table 1. Re-calculated properties of the compiled low- metallicity galaxies (continued) 



Object 


R23 


n H (S+) a 


t(0 2 +) b 


12 + log(0/H) 


Ref. c 


SBS 1129+576a 


3.851 


± 


0.111 


< 248 


1.899+°?" 


7.369 


± 


0.069 


G03a 


SBS 1129+576b 


5.804 


± 


0.247 


< 287 


2.094+°-?^ 

KJtJ — 0.275 


7.475 


± 


0.064 


G03a 


SBS 1159+545 


5.620 


± 


0.051 


57 +B * 

^ ' —50 


1.852 ± 0.020 


7.491 


± 


0.009 


198 


SBS 1205+557 


7.067 


± 


0.103 


< 76 


1.607 ± 0.067 


7.752 


± 


0.029 


197 


SBS 1211+540 


6.814 


± 


0.073 


168+! 24 


1.699 ± 0.024 


7.644 


± 


0.013 


194 


SBS 1222+614 


9.071 


± 


0.050 


r>r> + 28 

— 22 


1.425 ± 0.012 


7.951 


± 


0.009 


197 


SBS 1249+493 


7.359 


± 


0.087 


< 10 


1.648 ± 0.024 


7.721 


± 


0.012 


198 


SBS 1319+579 A 


9.913 


± 


0.066 


145+25 


1.310+°°" 

-0.012 


8.084 


± 


0.010 


197 


SBS 1319+579 B 


6.797 


± 


0.189 


40+" 3 

w — 40 


—0.175 


7.911 


± 


0.102 


197 


SBS 1319+579 C 


7.054 


± 


0.061 


20+S 

w — 20 


1.136+°°^ 

—0.036 


8.136 


± 


0.032 


197 


SBS 1331+493 


8.129 


± 


0.110 


164+«? 

— 81 


1.602 + 0.027 


7.780 


± 


0.016 


194 


SBS 1331+493S 


6.308 


± 


0.138 


< 79 


1 oro+0.086 
1,JJ,J -0.0S8 


7.885 


± 


0.056 


T95 


SBS 1415+437 


5.677 


± 


0.025 


65+^ 

yj ~' — 30 


1.703 ± 0.011 


7.586 


± 


0.005 


198 


SBS 1415+437el 


5.649 


± 


0.025 


48+i 1 , 

— 30 


1.657 + 0.010 


7.601 


± 


0.005 


G03c 


SBS 1415+437e2 


5.290 


± 


0.079 


79+" 4 

— 79 


1.597 ± 0.056 


7.614 


± 


0.028 


G03c 


SBS 1420+544 


9.683 


± 


0.070 


< 10 


1.764 + 0.011 


7.752 


± 


0.006 


198 


SBS 1533+469 


8.690 


± 


0.226 


46+ 4 ° 6 


l-383l°;°55 


7.984 


± 


0.035 


T95 


SBS 1533+574 A 


7.507 


± 


0.081 




1.444 ± 0.057 


7.883 


± 


0.031 


197 


SBS 1533+574 B 


9.107 


± 


0.083 


< 23 


1.246 ± 0.029 


8.124 


± 


0.023 


197 


SDSS J0113+0052 


3.402 


± 


0.211 


(10) d 


01 74-0.325 
— 301 


7.163 


± 


0.076 


106 


SDSS J0519+0007 


6.176 


± 


0.122 


q 7 O + 308 

O ( O 220 


2.078 ± 0.036 


7.420 


± 


0.015 


I04b 


SDSS J2104-0035 N 


4.044 


± 


0.093 


(10) d 


2.008 ± 0.066 


7.257 


± 


0.028 


106 


UGC 4305 5 


5.383 


± 


0.354 


(10) d 


1.607™ 

— U. 140 


7.647 


± 


0.064 


L03b 


UGC 4305 7 


5.150 


± 


0.280 


(10) d 


1.354+g- 2 2 ° 2 ! 


7.806 


± 


0.129 


L03b 


UGC 4305 8 


5.010 


± 


0.306 


(10) d 


1.514+°'!^ 

— U. 193 


7.677 


± 


0.086 


L03b 


UGC 4305 9 


5.371 


± 


0.379 


(10) d 


1.486+°"* 


7.699 


± 


0.072 


L03b 


UGC 4483 


4.795 


± 


0.052 


' — 72 


1.657 ± 0.026 


7.540 


± 


0.012 


194 


UGC 6456 


5.918 


± 


0.062 


9Q+52 

— 29 


1.547 ± 0.022 


7.696 


± 


0.012 


197 


UGC 6456 1 


4.395 


± 


0.399 


(10) d 


2.089+°!^ 

— 0.436 


7.355 


± 


0.108 


L03b 


UGC 6456 2 


5.180 


± 


0.372 


(10) d 


1.768+°J" 


7.519 


± 


0.077 


L03b 


UGC 9128 


4.178 


± 


0.166 


198 ± 16 


1.320+ ,'! 2 ,* 

0.127 


7.745 


± 


0.080 


L03b 


UGC 9497 c 


7.138 


± 


0.101 


< 77 


1.796 ± 0.030 


7.608 


± 


0.015 


G03b 


UGC 9497 e 


4.451 


± 


0.307 


< 453 


1.657+ ,' 43 * 

1 —0.440 


7.524 


± 


0.172 


G03b 


UM 133 


6.784 


± 


0.112 


< 26 


1.676 ± 0.032 


7.692 


± 


0.014 


I04b 


UM 238 


10.786 


± 


0.219 


288+K? 

— ol 


1.250 ± 0.016 


8.177 


± 


0.017 


I04b 


UM 311 


7.075 


± 


0.077 


75+ 4 ° 


0.977+°°^ 

yj.u i i —0.039 


8.374 


± 


0.044 


198 


UM 396 


9.391 


± 


0.177 


o 7 +44 
— 37 


1.136 + 0.016 


8.238 


± 


0.019 


I04b 


UM 420 


7.814 


± 


0.182 


< 79 


1 007+0. 081 
L.OOI _q 0g3 


7.941 


± 


0.049 


198 


UM 422 


10.519 


± 


0.200 


< 57 


1.296 + 0.014 


8.121 


± 


0.014 


I04b 


UM 439 


11.735 


± 


0.226 


1771" 


1.411 + 0.015 


8.069 


± 


0.013 


I04b 


UM 448 


6.225 


± 


0.067 


138+1 


1 99n +0.059 


8.018 


± 


0.047 


198 


UM 461 


8.518 


± 


0.200 




1.615 + 0.042 


7.782 


± 


0.027 


198 


UM 462 SW 


8.283 


± 


0.052 


< 10 


1.378 + 0.015 


7.960 


± 


0.010 


198 



a Gas density of the S+ regions in units of cm 3 . 

b Gas temperature of the 2 + regions in units of 10 4 K. 

c References. — GOO: Guseva et al. (2000), G03a: Guseva et al. (2003a), G03b: Guseva et al. (2003b), G03c: Guseva et al. 
(2003c), 194: Izotov et al. (1994), 197: Izotov et al. (1997), 198: Izotov & Thuan (1998), I04a: Izotov et al. (2004), I04b: 
Izotov & Thuan (2004), 106: Izotov et al. (2006a), K96: Kobulnicky & Skillman (1996), L03a: Lee et al. (2003a), L03b: Lee 
et al. (2003b), L04: Lee et al. (2004), M04: Melbourne et al. (2004), P92: Pagel et al. (1992), P06: Papaderos et al. (2006), 
T95: Thuan et al. (1995), V03: Vflchez & Iglesias-Paramo (2003). 

d Flux ratio of [Sn] is not given in literature. We adopt nn(S+) = 10 cm~ 3 to calculate the oxygen abundance for these 
objects. 
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Sample 




Median of 12+log(0/H) 


Mean of 12+log(0/H) 


RMS of 12+log(0/H) 


Sample A 


139 


8.010 


8.003 


0.233 


Sample B 


120 


7.936 


7.858 


0.303 


Sample C 


48497 


9.016 


8.976 


0.166 


Sample A+B 


259 


7.980 


7.936 


0.277 



Table 3. Means and RMSs of emission-line flux ratios of the galaxies in the sample A+B a 



Oxygen Abundance 



logife log«|p 



i F([OIH]A5007) 
1U 6 F([NII]A6584) 



log ££NnA6584]_ 



F( Oil A3727) 



, F([NII]A6584) 
1U & F([SII]A6720) 



i F([OIII]A5007) 
1U & F([OII]A3727) 



7.05<12+log(O/H)<7.45 
7.45 <12+log(0/H)< 7.55 
7.55 <12+log(0/H)< 7.65 
7.65 <12+log(0/H)< 7.75 
7.75 <12+log(0/H)< 7.85 
7.85 <12+log(0/H)< 7.95 
7.95 <12+log(0/H)< 8.05 
8.05 <12+log(0/H)< 8.15 
8.15 <12+log(0/H)< 8.25 
8.25 <12+log(0/H)< 8.35 
8.35 <12+log(0/H)< 8.55 



0.601 
(0.108) 

0.780 
(0.084) 

0.809 
(0.068) 

0.843 
(0.039) 

0.909 
(0.039) 

0.928 
(0.045) 

0.941 
(0.042) 

0.944 
(0.053) 

0.938 
(0.044) 

0.917 
(0.036) 

0.892 
(0.034) 



-2.452 
(0.238) 

-2.122 
(0.272) 

-1.844 
(0.195) 

-1.846 
(0.247) 

-1.887 
(0.205) 

-1.624 
(0.197) 

-1.578 
(0.213) 

-1.481 
(0.177) 

-1.375 
(0.142) 

-1.415 
(0.196) 

-1.124 
(0.176) 



2.415 
(0.196) 

2.278 
(0.343) 

2.006 
(0.274) 

2.050 
(0.290) 

2.162 
(0.260) 

1.884 
(0.262) 

1.858 
(0.280) 

1.759 
(0.255) 

1.631 
(0.185) 

1.641 
(0.260) 

1.259 
(0.221) 



-1.657 
(0.269) 

-1.477 
(0.120) 

-1.376 
(0.183) 

-1.431 
(0.176) 

-1.458 
(0.129) 

-1.281 
(0.167) 

-1.262 
(0.179) 

-1.227 
(0.149) 

-1.119 
(0.114) 

-1.201 
(0.101) 

-1.084 
(0.109) 



-0.668 
(0.178) 

-0.606 
(0.111) 

-0.566 
(0.197) 

-0.599 
(0.075) 

-0.621 
(0.112) 

-0.531 
(0.140) 

-0.484 
(0.118) 

-0.450 
(0.109) 

-0.348 
(0.120) 

-0.333 
(0.064) 

-0.173 
(0.122) 



0.920 
(0.247) 

0.814 
(0.262) 

0.652 
(0.356) 

0.492 
(0.233) 

0.719 
(0.233) 

0.579 
(0.307) 

0.571 
(0.218) 

0.511 
(0.247) 

0.514 
(0.168) 

0.436 
(0.288) 

0.258 
(0.087) 



Mean and RMS of each emission-line flux ratio are given in the upper and lower rows. RMSs are given in parenthesis. 
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Table 4. Means and RMSs of emission-line flux ratios of the galaxies in sample C a 



Oxygen Abundance 



logi? 2 



, F([NII]A6584) 

10 & F(H«) 



! F([OIII]A5007) 

10 6 F([NII]A6584) 



logi 



F( 


Nil 


A6584) 


F( 


Oil 


A3727) 



I , F( MI Al,r,84) 
10 6 F([SII]A6720) 



, F^[OIU]A50072 
10 & F([OII]A3727) 



8.15 <12+log(0/H)< 8.25 
8.25 <12+log(0/H)< 8.35 
8.35 <12+log(0/H)< 8.45 
8.45 <12+log(0/H)< 8.55 
8.55 <12+log(0/H)< 8.65 
8.65 <12+log(0/H)< 8.75 
8.75 <12+log(0/H)< 8.85 
8.85 <12+log(0/H)< 8.95 
8.95 <12+log(0/H)< 9.05 
9.05 <12+log(0/H)< 9.15 
9.15 <12+log(0/H)< 9.25 



0.835 
(0.043) 

0.825 
(0.053) 

0.817 
(0.060) 

0.805 
(0.061) 

0.772 
(0.052) 

0.711 
(0.052) 

0.637 
(0.052) 

0.559 
(0.058) 

0.454 
(0.071) 

0.324 
(0.084) 

0.168 
(0.074) 



-1.138 
(0.081) 

-1.031 
(0.077) 

-0.934 
(0.088) 

-0.851 
(0.119) 

-0.812 
(0.101) 

-0.689 
(0.083) 

-0.596 
(0.068) 

-0.516 
(0.058) 

-0.447 
(0.056) 

-0.415 
(0.052) 

-0.398 
(0.056) 



1.102 
(0.171) 

0.936 
(0.187) 

0.768 
(0.199) 

0.644 
(0.237) 

0.572 
(0.216) 

0.300 
(0.188) 

0.060 
(0.160) 

-0.159 
(0.139) 

-0.383 
(0.123) 

-0.591 
(0.126) 

-0.706 
(0.192) 



-1.196 
(0.077) 

-1.113 
(0.082) 

-1.056 
(0.084) 

-0.985 
(0.102) 

-0.910 
(0.091) 

-0.780 
(0.084) 

-0.642 
(0.075) 

-0.508 
(0.077) 

-0.345 
(0.088) 

-0.174 
(0.100) 

0.016 
(0.085) 



-0.532 
(0.042) 

-0.480 
(0.038) 

-0.444 
(0.033) 

-0.397 
(0.034) 

-0.334 
(0.047) 

-0.238 
(0.051) 

-0.146 
(0.051) 

-0.053 
(0.052) 

0.054 
(0.057) 

0.171 
(0.057) 

0.286 
(0.044) 



-0.089 
(0.140) 

-0.191 
(0.170) 

-0.283 
(0.174) 

-0.309 
(0.184) 

-0.329 
(0.170) 

-0.466 
(0.147) 

-0.573 
(0.135) 

-0.663 
(0.120) 

-0.731 
(0.108) 

-0.777 
(0.101) 

-0.758 
(0.118) 



Mean and RMS of each emission-line flux ratio are given in the upper and lower rows. RMSs are given in parenthesis. 



Table 5. Coefficients of the best- fit polynomials for the observed relations between the emission- line flux ratios and 
the oxygen abundance, where logi? = a + a\x + a 2 x 2 + a 3 x 3 [x = \og{Z / Zq) = 12+log(0/H)-8.69]. 



Flux ratio (logi?) 


ao 


ai 


0,2 


as 


log i?23 


+7.1806E-1 


-6.9548E-1 


-6.2220E-1 


-6.3169E-2 


log [F([Nii]A6584)/F(Hq)] 


-6.8307E-1 


+8.9881E-1 


-5.2302E-1 


-2.2040E-1 


log [F( [Om] A5007) /F( [Nil] A6584)] 


+3.2921E-1 


-2.2578E+0 


-4.1699E-2 


+3.7941E-1 


log [F( [Nil] A6584) /F( [On] A3727)] 


-7.9322E-1 


+1.1399E+0 


+7.8929E-1 


+2.7101E-1 


log [F( [Nn] A6584) /F( [Sn] A6720)] 


-2.5214E-1 


+7.4100E-1 


+5.8181E-1 


+1.7963E-1 


log [F([Oiii]A5007)/F([Oii]A3727)] 


-3.0777E-1 


-1.1210E+0 


-1.4359E-1 





Nagao et al.: Gas metallicity diagnostics 



27 



Table 6. Coefficients of the best- fit polynomials for the observed relations between the emission-line flux ratios and 
the oxygen abundance, where logi? = &o + biy + b 2 y 2 + b 3 y 3 [y = 12+log(0/H)]. 



Flux ratio (logJ?) 



60 



61 



b 2 



63 



log R 23 +1.2299E+0 -4.1926E+0 +1.0246E+0 -6.3169E-2 

log [F([Nn]A6584)/F(Ha)] +9.6641E+1 -3.9941E+1 +5.2227E+0 -2.2040E-1 

log [F([Oiii]A5007)/F([Nii]A6584)] -2.3218E+2 +8.4423E+1 -9.9330E+0 +3.7941E-1 

log [F([Nii]A6584)/F([Oii]A3727)] -1.2894E+2 +4.8818E+1 -6.2759E+0 +2.7101E-1 

log [F([Nn]A6584)/F([Sn]A6720)] -8.0632E+1 +3.1323E+1 -4.1010E+0 +1.7963E-1 

log [F([Oiii]A5007)/F([Oii]A3727)] -1.4089E+0 +1.3745E+0 -1.4359E-1 



Table 7. Means and RMSs of additional emission-line flux ratios of the galaxies in the sample A+B a 



Oxygon Abundance 


]nr ,f([NeIIIlA3869) 
s F([OII]A3727) 


, F(Ha+[NII]AA6548,6584) 
S F([SII]A6720) 


l nr ^([Oin]A5()()7) 
1U & F(H/3) 


W F([OIIlA3727) 
1U & F(H/3) 


7 05 < 1 9-t-liWO /HI ^74^ 


-0.092 


1.737 


0.421 


-0.486 




CO 144) 


(0.156) 


(0.134) 


(0.186) 


7.45 <12+log(0/H)< 7.55 


-0.285 


1.481 


0.594 


-0.210 




(0.255) 


(0.209) 


(0.100) 


(0.161) 


7.55 <12+log(0/H)< 7.65 


-0.426 


1.254 


0.600 


-0.050 




(0.368) 


(0.256) 


(0.124) 


(0.226) 


7.65 <12+log(0/H)< 7.75 


-0.545 


1.188 


0.653 


0.127 




(0.202) 


(0.175) 


(0.062) 


(0.193) 


7.75 <12+log(0/H)< 7.85 


-0.374 


1.217 


0.724 


0.003 




(0.219) 


(0.177) 


(0.072) 


(0.171) 


7.85 <12+log(0/H)< 7.95 


-0.504 


1.038 


0.708 


0.140 




(0.287) 


(0.171) 


(0.082) 


(0.216) 


7.95 <12+log(0/H)< 8.05 


-0.510 


1.036 


0.731 


0.163 




(0.222) 


(0.149) 


(0.082) 


(0.144) 


8.05<12+log(O/H)<8.15 


-0.592 


0.992 


0.732 


0.216 




(0.228) 


(0.150) 


(0.092) 


(0.163) 


8.15<12+log(0/H)<8.25 


-0.603 


0.968 


0.712 


0.215 




(0.150) 


(0.101) 


(0.069) 


(0.103) 


8.25 <12+log(0/H)< 8.35 


-0.694 


1.038 


0.684 


0.244 




(0.346) 


(0.197) 


(0.086) 


(0.188) 


8.35 <12+log(0/H)< 8.55 


-0.900 


0.932 


0.590 


0.361 




(0.095) 


(0.082) 


(0.060) 


(0.074) 



a Mean and RMS of each emission-line flux ratio are given in the upper and lower rows. RMSs are given in parenthesis. 
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Table 8. Means and RMSs of additional emission-line flux ratios of the galaxies in sample O 



Oyvcfpn A lm nH a n pp 




l nr ,f([NeIII]A3869) 
1U & F([OII]A3727) 


,„„F(Ha+[NII]AA654S,6584) 
10 S F([SII]A6720) 


, F([OIII]A5007) 
1U & F(H/3) 


, F([OII]A3727) 
1U & F(H/3) 


8.15<12+log(0/H)<8 


.25 


-1.060 


0.637 


0.405 


0.512 






(0.112) 


(0.081) 


(0.102) 


(0.062) 


8.25<12+log(0/H)<8 


35 


-1.084 


0.591 


0.357 


0.549 






(0.117) 


(0.089) 


(0.119) 


(0.079) 


8.35<12+log(0/H)<8 


45 


-1.132 


0.540 


0.281 


0.580 






(0.111) 


(0.098) 


(0.123) 


(0.091) 


8.45<12+log(0/H)<8 


55 


-1.113 


0.518 


0.250 


0.578 






(0.109) 


(0.105) 


(0.134) 


(0.086) 


8.55<12+log(0/H)<8 


65 


-1.190 


0.545 


0.219 


0.555 






(0.121) 


(0.087) 


(0.129) 


(0.075) 


8.65<12+log(0/H)<8 


75 


-1.301 


0.538 


0.070 


0.544 






(0.121) 


(0.068) 


(0.124) 


(0.062) 


8.75<12+log(0/H)<8 


85 


-1.339 


0.553 


-0.077 


0.501 






(0.121) 


(0.057) 


(0.120) 


(0.057) 


8.85<12+log(0/H)<8 


95 


-1.400 


0.583 


-0.218 


0.446 






(0.123) 


(0.048) 


(0.116) 


(0.060) 


8.95<12+log(0/H)<9 


05 


-1.371 


0.639 


-0.369 


0.352 






(0.266) 


(0.043) 


(0.116) 


(0.074) 


9.05<12+log(O/H)<9 


15 




0.727 


-0.530 


0.228 






(-) 


(0.047) 


(0.119) 


(0.089) 


9.15<12+log(0/H)<9 


25 




0.823 


-0.611 


0.078 






(-) 


(0.043) 


(0.177) 


(0.085) 



Mean and RMS of each emission-line flux ratio are given in the upper and lower rows. RMSs are given in parenthesis. 



Table 9. Coefficients of the best-fit polynomials for the observed relations between the additional emission- line flux 
ratios and the oxygen abundance, where logi? = a + a\x + a 2 x 2 + a 3 x 3 [x = \og(Z/Z Q ) = 12+log(0/H)-8.69]. 



Flux ratio (logi?) 




ai 


a 2 


0-3 


log [F( [Nem] A3869)/F( [On] A3727)] 


-1.2547E+0 


-7.0929E-1 


+3.0497E-1 


+1.6784E-1 


log [F(Ha+[Nii]AA6548,6584)/F([Sn]A6720)] 


+5.6097E-1 


-7.9971E-2 


+9.8562E-1 


+3.4069E-1 


log [F([Om]A5007)/F(H/3)] 


+1.6366E-1 


-1.3785E+0 


-8.4778E-1 


+9.1853E-3 


log [F([On]A3727)/F(H/3)] 


+5.3481E-1 


-2.0792E-1 


-1.1353E+0 


-3.5951E-1 



Table 10. Coefficients of the best-fit polynomials for the observed relations between the additional emission-line flux 
ratios and the oxygen abundance, where logi? = bo + b\y + b2y 2 + b^y 3 [y = 12+log(0/H)]. 



Flux ratio (logi?) 


bo 


bi 


b 2 


63 


log [FQNem] A3869)/F([On] A3727)] 


-8.2202E+1 


+3.2014E+1 


-4.0706E+0 


+1.6784E-1 


log [F(Ha+[Nn]AA6548,6584)/F([Sn]A6720)] 


-1.4789E+2 


+5.9974E+1 


-7.8963E+0 


+3.4069E-1 


log [F([Om]A5007)/F(H/3)] 


-5.7906E+1 


+1.5437E+1 


-1.0872E+0 


+9.1853E-3 


log [F([On]A3727)/F(H/3)] 


+1.5253E+2 


-6.1922E+1 


+8.2370E+0 


-3.5951E-1 



